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ABSTRACT 


The  rapidly  expanding  fields  of  digital  control  and  advanced  power  electronic  circuitry 
enable  researchers  and  designers  to  investigate  novel  approaches  of  providing  shipboard 
service  power  requirements.  The  Shipboard  More-Electric  concept  is  based  on  using  high 
efficiency  DC-DC  PWM  converter  modules  in  place  of  existing  low  efficiency  high 
maintenance  power  distribution  systems  in  the  future  generation  Navy  combatant  design. 

This  thesis  describes  the  design  and  implementation  of  a  reduced-order  model  of  a 
Primary  Ship  Service  Converter  Module  (PSSCM).  Using  the  Advanced  Continuous 
Simulation  Language  (ACSL),  the  PSSCM  is  simulated  to  predict  performance  and 
stability.  Various  voltage  sources,  loads  and  controllers  are  studied  to  ensure  the  suitability 
and  reliability  of  the  system. 
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L  INTRODUCTION 

As  technology  advances,  new  and  better  methods  are  being  researched  and  developed 
for  providing  higher  performance  and  realizing  cost  savings  in  shipboard  power  systems. 
In  an  attempt  to  improve  reliability,  enhance  survivability,  and  broaden  both  commonality 
and  flexibility,  conventional  circuit  topologies  are  being  replaced  by  state  of  the  art  high 
efficiency  PWM  converters. 

The  latest  concept  of  distributing  DC  pow^  for  both  the  ship’s  sovice  power  and  the 
ship’s  prc^ulsion  systems  is  known  as  the  Integrated  Power  System  (IPS).  The  DC  powo* 
is  delivered  to  all  zones  in  the  ship  through  a  port  and  starboard  bus  from  the  main  turbine 
generators.  Through  isolating  the  various  zones  so  that  there  are  no  interzone  electrical 
connections,  the  generic  description  Zonal  Electric  Distribution  (ZED)  is  arrived  at.  The 
Primary  Ship  Service  Converte*  Mcxlule  (PSSCM)  is  the  key  to  providing  eacdi  zone  with 
regulated  and  stable  power.  Each  zone  contains  both  DC  and  AC  load  requirements.  A 
representation  of  a  PSSCM,  upon  employing  various  assumptions  and  reductions,  can  be 
simulated  using  the  Advanced  Continuous  Simulation  Language  (ACSL).  The  PSSCM  in 
the  reduced  form  is  a  DC-DC  Buck  converter.  The  concepts  and  operations  of  the  Buck 
converter  are  analyzed  and  applied  to  simulate  the  PSSCM  which  is  then  integrated 
together  with  some  typical  zone  loads.  Several  studies  are  analyzed  to  determine  the 
stability  and  performance  of  the  PSSCM.  Fixed  constant  sources  and  active  sources  are 
specified  to  investigate  and  critique  the  perfmmance  of  the  PSSCM. 

To  enhance  the  stability  of  the  PSSCM,  several  controllers  are  implemented  to  control 
the  output  voltage  through  modulation  of  the  duty  cycle.  From  the  three  voltage  controllers 
analyzed,  it  is  shown  tiiat  the  multi-loop  controller  had  the  most  desirable  performance.  The 
controllM"  increases  the  quality  of  the  output  voltage  by  reducing  the  oscillations,  voltage 
spikes  and  response  times  during  transients. 

In  order  to  assess  zone  power  and  transient  recovery  requirements  fw  SSCM’s,  a 
number  of  simulations  are  performed  including  changes  in  load  resistance  and  the  start-up 
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of  an  inverter-driven  induction  machine.  Simulations  are  also  perfcamed  to  study  the 
effects  of  two  PSSCMs  running  simultaneously  from  the  same  poww  source,  effectively 
modeling  a  two  zone  configuration.  A  normal  destroyer-size  combatant  may  incorporate 
eight  or  more  electrical  zones.  The  final  issue  considoed  is  a  paralleled  combination  of 
Buck  converters  where  the  control  signals  are  interlaced  in  an  attempt  to  exploit  further 
control  of  the  ou^ut  voltage. 
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n.  THE  MORE-ELECTRIC  SYSTEM 

A.  THE  MORE-ELECTRIC  CONCEPT 

The  Navy’s  move  toward  more-electric  is  an  attempt  to  find  alternatives  to  the  present 
distribution  architectures  which  provide  ship  propulsion  and  ship  service  power 
requirements.  The  integrated  nature  of  this  endeavor  stems  from  the  reqiiirement  that  the 
ship’s  sowice  power  must  be  derived  from  the  same  source  as  the  ship’s  main  propulsion. 
It  is  believed  that  this  tq>ology  will  provide  ina*eased  efficiency  and  cost  reductions.  By 
minimizing  the  number  of  piece  parts  and  by  using  pre-designed  modules  with  standard 
components,  tiie  More-Electric  ship  alternative  is  intended  to  reduce  the  acquisition  and  life 
(^cle  costs.  In  addition,  it  is  the  goal  of  the  more-electiic  ship  to  maintain  the  poformance 
c^abilities  of  combat  vessels  and  readily  facilitate  future  modifications  and  upgrades[l]. 

B.  THE  MORE-ELECTRIC  ARCHITICTURE 

1.  Electric  Propulsion 

The  generate  is  driven  by  a  22  MW  Intercooled  Recuperative  (ICR)  gas  turbine 
engine  at  3600  ipm.  This  is  the  next  generation  gas  turbine  which  the  Navy  is  developing 
to  replace  ot  upgrade  the  LM2S00.  The  ICR  turbine  offos  superior  transient  performance 
over  the  LM2500  without  an  increase  in  size  ca*  weight.  The  turbine  drives  a  permanent 
magnet  synchronous  generator  with  multiple  phase  windings  Each  winding  has  an 
independent  rectifier  circuit  which  supplies  either  the  DC  ship  propulsion  or  the  DC  ship 
service  bus. 

2.  Direct  Current  2^nal  Electric  Distribution  (DC  ZED) 

Zonal  Electric  Distribution  (ZED)  is  a  system-oriented  approach  for  supplying 
power  to  the  primary  and  secondary  loads  present  on  a  typical  Navy  ship.  The  DC  7Tn~)  was 
selected  over  an  AC  ZED  because  of  the  available  split-plant  configuration,  the  near 
instantaneous  power  transfer,  the  added  control  flexibility,  and  the  inherent  current-limit 
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protection  built  into  the  power  electronic  equipment.  Due  to  advancements  in  poww 
electronics  which  include  greater  device  current  densities,  increased  reliability  and 
efficiency,  and  reduction  in  costs,  DC  distribution  is  now  a  feasible  alternative  to 
conventional  AC  systems[2].  The  disadvantages  inherent  to  DC  ZED  are  arcing,  fault 

detection  (no  zero  crossing),  and  isolation. 

The  DC  is  sectioned  into  a  port  and  a  starboard  bus  with  each  bus  linking  a 
zone  through  a  Primary  Shipboard  Service  Converter  Module  (PSSCM).  The  design  of  two 
buses,  one  on  the  port  and  the  other  on  the  starboard,  provides  the  ship  with  redundancy  in 
case  of  repair  or  ship  casualty.  The  PSSCM’s  supply  the  appropriate  DC  voltage  to  the  Ship 
Service  Inverter  Module  (SSIM’s)  and  Ship  Service  Converter  Module  (SSCM’s)  in  each 
zone.  The  SSIM’s  deliver  power  to  the  AC  loads,  while  the  SSCM’s  supply  the  DC  loads 
in  each  zone.  Figure  2-1  illustrates  flie  architecture  of  the  More-Electric  System  [1]. 

C.  RESEARCH  STATUS 

The  Naval  Surface  Warfare  Center  (NSWC)  at  Annqrolis,  Maryland  and  the  Power 
Systems  Group(PSG)  of  Power  Paragon  Inc.  have  combined  efforts  to  produce  a  Primary 
Ship  Service  Converter  Module.  The  entire  project,  including  design,  manufacturing  and 
the  SSCM  will  be  a  joint  effwt.  The  PSSCM  is  a  DC-DC  Buck  type  converter. 

This  PSSCM  win  be  designed  and  a  prototype  developed,  installed,  and  operated  in  the 
laboratory  of  the  NSWC  fw  DC  ZED's  system  performance  evaluations.  PSG  is  presently 
working  on  the  power  pcxtion  design  of  the  PSSCM  system.  NSWC  and  NPS  are  presently 
researching  the  control  portion  of  the  PSSCM  system.  The  Advanced  Surface  Machinery 
Programs  (ASMP)  office  at  NSWC  has  designed  a  Reduced  Scale  Advanced  Development 
(RSAD)  model  fw  the  purpose  of  promoting  the  DC  ZED.  A  Full  Scale  Advanced 
Development  (FSAD)  model  will  facilitate  the  validation  and  appropriate  modification  of 
the  computer  models  presently  being  used  to  assess  stability  and  performance 
characteristics  [1],[3]. 
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Figure  2.1:  The  More-Elechic  System. 
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in.  DC-DC  BUCK  CONVERTER 

A.  WHY  THE  DC  BUCK  CONVERTER 

The  biggest  consideration  in  the  use  of  step-down  converters  over  conventional  linear 
control  circuits  is  efficiency.  Converters  can  operate  at  efficiencies  greater  than  90%,  while 
conventional  voltage  divider  circuits  operate  with  mudi  largo*  power  losses,  less  efficiency 
and  have  no  capacity  to  regulate  the  output  voltage. 

The  four  basic  Switched-Mode  PWM  DC-DC  convoters  are  the  Buck,  Boost,  Buck- 
Boost  and  Cuk.  The  Buck  converter  steps  down  a  large  voltage  to  a  usable  lower  level.  The 
Boost  converter  as  the  name  implies,  ino-eases  an  existing  voltage,  while  the  Buck-Boost 
and  Cuk  can  eitho*  step-up  or  step-down  a  given  source  to  a  desired  voltage  level.  The  Cuk 
converter  is  an  inverse  dual  of  the  Buck-Boost  converter. 

In  the  M(X'e-Electric  shipboard  design,  the  voltage  level  will  be  at  a  mayimntn  level 
when  it  leaves  the  integrated  generatorAectiiler  unit.  After  being  distributed,  the  DC 
voltage  must  then  be  reduced  within  each  zone  by  a  PSSCM.  In  each  zone,  the  voltage  may 
be  further  reduced  by  additional  SSCM’s.  It  is  apparent  that  the  More-Electric  shipboard 
design  requires  only  step-down  converters,  therefore  the  Buck  converter  is  the  best  choice 
for  the  PSSCM. 

B.  GENERAL  BUCK  CONVERTER  CONFIGURATION  AND  OPERATION 

The  basic  DC  Buck  converter  is  shown  in  Figure  3-1.  El  is  the  input  voltage.  The 

ouq)ut  voltage  across  the  capacitor  is  Vcl.  The  capacitor  is  assumed  to  be  large  so  the 
AC  portion  of  the  voltage  is  relatively  small  and  can  be  neglected  for  various  purposes  of 
analysis.  The  diode,  Dl,  is  assumed  to  be  an  ideal  diode  and  thus  has  zero  voltage  drop 
when  forward  biased  and  carries  zero  current  when  reverse  biased.  The  usual  0.7  to  2.0  volt 
diode  voltage  drop  can  be  reasonably  neglected  in  the  remainder  of  the  analysis  because  the 
circuit  voltages  imder  consideration  are  much  larger.  The  placement  of  the  diode  in  the 
circuit  provides  a  path  for  the  inductor  current  when  the  switch  is  open.  The  switch,  SI,  is 
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assumed  to  be  an  ideal  switch  with  two  states,  open  or  closed.  If  the  switch  is  open  and  D1 
is  conducting,  the  voltage  awoss  the  switch  is  El  volts.  When  the  switch  is  closed,  it  is 
assumed  that  the  voltage  aaoss  the  switch  instantaneously  transitions  from  El  volts  to  zero 
volts.  Switching  transients  and  switch  power  loss  are  ignored  for  convenience.  The  switch 
is  turned  on  and  off  at  a  frequency  generally  between  5  kHz  and  20  kHz  for  high-power 
£q)plications.  The  switching  process  introduces  a  ripple  in  the  ouqjut  voltage.  For  the 
analysis  and  modeling  presented  in  this  work,  harmonic  losses  and  power  factor 
considerations  are  assumed  to  be  negligible. 

Even  though  the  switching  frequency  is  constant,  the  ratio  between  the  switch  on  and 
off  time  will  vary  according  to  some  prescribed  control  law.  The  ratio  of  the  switch  on-time 
to  the  switching  period  is  termed  the  duty  cycle  without  any  other  consideration,  this 
p^odic  switching  will  introduce  considerable  harmonics  into  the  rest  of  the  circuit.  The 
positioning  of  the  Buck  chopper  inductor  and  ciq)acitor  establishes  a  low-pass  filter  with 
cutoff  frequency  Wo=  l/sqrt(Ll*Cl)  which  filters  out  much  of  the  switching  harmonics 
before  they  reach  the  ouq)ut  load,  Rl. 


Figure  3.1;  DC  Buck  Configuration 
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1.  Switch-Closed  Configuration 

When  the  switch  is  closed  the  Buck  Chopper  circuit  topology  is  given  by  Figure 
3.2.  In  this  mode,  the  diode  is  reverse  biased  and  can  be  ignored. 

By  Kirchoff  s  voltage  law,  the  following  equations  are  derived: 


E1=Vli  +Vci 

(3.1) 

El=Ll*p(iLl)+Vci 

(3.2) 

p(iLi)  =  (El-Vci)/Ll 

(3.3) 

Imax-Imin  =  ((El-Vci)/L1)*(D*T) 

(3.4) 

Note  that  p  is  the  Heaviside  operate,  T  is  the  switch  poiod  and  equal  to  the 
inverse  frequency  and  D,  the  duty  cycle,  is  the  portion  of  the  time  the  switch  is  on  during 
each  period  (0<D<1). 


Figure  3.2:  Switch-Qosed  Configuration 


9 


2.  Switch-Open  Configuration 

When  the  switch,  SI,  is  open  the  Buck  Chopper  circuit  t(^ology  is  given  by 
Figure  3.3. 

Again  by  Kirchoffs  voltage  law,  the  following  equations  are  derived 


0=Vli+Vci 

(3.5) 

0=Ll*p(iLi)  +Vci 

(3.6) 

p(iLl)=  ~^Clf  LI 

(3.7) 

Imax-hnin  =(V  qi/  L1)*(D’*T) 

(3.8) 

Note  that  D’  is  the  complement  of  the  duty  cycle  and  represents  the  ratio  of  time 
that  the  switch  is  off  to  the  switching  freqency  (D  +  D’  =1).  If  (3.4)  and  (3.8)  are 
simultaneously  solved,  the  following  equation  is  derived: 

Vci=D*El  (3.9) 

Note  that  the  ouq>ut  voltage,  Vcl,  is  proportional  to  the  input  voltage,  El,  and  is 
directly  adjusted  by  the  duty  cycle.  Equation  (3.9)  is  valid  if  all  components  are  ideal,  wire 
loss  is  negligible,  inductive  resistance  is  neglected,  and  continuous  inductor  current  is 
assumed. 


LI 


R1 
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C.  MODES  OF  OPERATION 


Li  the  previous  discussion  and  equation  derivations,  the  continuous  current  mode  of 
operation  was  assumed.  As  the  name  implies,  the  continuous  current  mode  occurs  when  the 
current  through  the  inductor  is  always  positive.  If  the  current  through  the  inductor  becomes 
zero  for  a  portion  of  the  switching  period,  then  the  converter  is  said  to  be  operating  in  the 
discontinuous  current  mode.  In  the  discontinuous  current  mode  the  converter  is  found  to  be 
less  stable,  and  higher  peak  currents  are  required  to  supply  the  same  average  ouq>ut  current 
as  in  the  continuous  current  mode.  As  a  result,  device  ratings  must  be  ino-eased  if  operation 
in  the  discontinuous  mode  is  desired. 

1.  Determining  the  Mode  of  Operation 

The  mode  of  operation,  whether  continuous  or  discontinuous,  depends  on  the 
switching  frequency,  the  load  resistance,  the  duty  cycle  and  the  inductor  size.  If  the  load 
resistance,  switching  frequency  and  duty  cycle  are  known  then  the  aitical  value  of 
inductance  can  be  found.  The  critical  inductor  value  is  the  determining  factor  in  assessing 
the  steady-state  operating  mode. 

For  continuous  conduction,  the  average  inductor  current  is  given  by 
iLx=(Iiiiax  +Imin)  /2  (3.10) 

The  current  through  the  inducts  is  equal  to  the  combined  cmrent  drawn  by  the 
cf^acitor  and  resists.  The  average  capacitor  current  is  zero  for  an  ideal  component  and  the 
current  through  the  load  is  equal  to  the  voltage  across  the  capacitor  divided  by  the  load 
resistance.  From  the  above,  it  follows  that 

iLl=k:i+iRl  (311) 

iLi=^  Vci/  R1  (3. 12) 

ff  (3.10)  is  substituted  in  for  iLi  in  (3.12),  and  (3.8)  is  used  to  substitute  for  Tmay 
while  Imin  is  set  to  zero,  the  following  expression  for  critical  inductance  can  be  derived: 

Ll>(D’*Rl*T)/2  (3.13) 
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Since  D‘  is  less  than  one,  LI  can  be  selected  greater  than  (Rl*T)/2  to  insure 
continuous  current  mode  operation. 

D.  PSSCM  SYSTEM 

The  primary  ship’s  service  convCTter  powo"  design,  illustrated  in  Figure  3.4,  can  be 
reduced  to  the  basic  Buck  converter  configuration  already  introduced  in  Figure  3.1.  The 
100  kw  PSSCM  was  designed  by  the  Power  System  Group.  Several  assumptions  are 
required  to  obtain  the  simplified  representation.  This  section  will  discuss  the  assumptions 
and  the  basic  component  design  of  the  PSSCM.  Though  there  are  several 
methodologies  that  may  be  employed  to  achieve  the  desired  converter  specifications,  the 
cmwpnnftnt  design  presented  here  is  based  on  the  Power  System  Group  analysis  [3]. 
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Figure  3.4:  PSSCM  Schematic  [3] 
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The  choice  of  inductor  size  for  the  Buck  converter  is  based  on  (3.14)  and  (3.15).  The 
change  in  current  is  equal  to  the  average  current  divided  by  the  time.  Recall,  time  is  the 
inverse  of  fiequency,  and  the  switching  frequency  was  selected  to  be  5  kHz.  Given  an 
ouq)ut  powra-  of  100  kW  and  an  output  voltage  of  750  volts,  the  rated  average  current  is 
found  to  be  133  amps.  The  values  of  the  minimum  and  maximum  currents  referred  to  in 
(3.10)  are  discussed  in  Chapter  V. 

nd  )  =  ?y  z=  — 133  _  565  000amp5/sec  (3.14) 

T  200e-6 

The  specifications  for  the  input  voltage  is  a  nominal  850  volts  with  an  admissible 
variation  of  +/-25  volts.  The  specification  for  the  output  voltage  is  750  volts. 
Conservatively,  the  maximum  voltage  anticipated  across  the  inductor  is  900  volts.  This  is 
a  conservative  value  since  the  maximum  input  voltage  is  specified  at  875  volts.  The 
maximum  voltage  across  the  inductor  occurs  at  start-up  when  the  mitial  condition  on  the 
capacitor  voltage  is  zero.  Solving  (3.3)  fw  LI  upon  substituting  for  the  maximxim  inductor 
voltage  yields 


LI 


Pdu) 


900 

665,000 


1.35mH 


(3.15) 


The  actual  drcuit  realization  uses  two  675uH  inductors  in  series,  instead  of  one 
1.35mH  inductor,  to  reduce  the  inductor  internal  capacitance  effects.  The  inductcxs  are 


identified  by  L2  and  L3  in  Figure  3.4. 

The  magnetic  energy  stored  in  the  inductor  is  denoted  as  eLi-  This  energy  is  a  function 
of  the  inductor  size  and  the  inductor  current.  In  particular  for  rated  current,  the  energy  is 
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The  standard  rule-of-thumb  specifies  that  the  minimum  capacitor  storage  energy 
should  be  on  the  ordo*  of  10  times  the  rated  inductor  storage  energy.  The  minimum 
cq>acit(»:  value  is  then  obtained  from 


Cmin  = 


2*eci  _  2»(10*12) 

2  2 

V  750 

c 


=  426uf 


(3.17) 


The  PSG  design  uses  a  capacitance  of  2600uf  which  is  a  factor  of  six  times  the 
minimum  required.  The  2600uf  capacitance  is  realized  by  combining  two  parallel 
combinations  of  three  3900  uf  ci^acitcvs  in  series.  These  capacitors  are  illustrated  in  Figure 
3.4  and  are  identified  by  C29  through  C31  and  C35  through  C37.  For  a  2600uf  output 
capacitor,  the  anticipated  peak-to-peak  ripple  voltage  is  0.1257  volts,  as  specified  by 

AW  _  (Imax-Imin)T 

- gJci 


The  switch,  SI,  is  an  Insulated  Gate  Bipolar  Transistor  (IGBT).  The  400  amp  1700 
volt  IGBT  is  gated  on  by  an  appropriate  driver,  not  illustrated  in  the  figure.  The  input  to  the 
gate  driver  establishes  the  duty  cycle  and  the  frequency  of  the  switch.  The  duty  cycle  is 
determined  by  the  type  of  controller  used.  A  through  discussion  of  controUo^  and 
controller  design  is  presented  in  Chq>te'  FV. 

The  main  Buck  converter  power  diode  is  identified  in  Figure  3.4  as  CRl.  To  reduce 
transient  voltage  ringing,  turn-off  switching  losses  and  IGBT  voltage  stresses,  two  snubber 
drcuits  are  implemented  across  the  IGBT.  For  identical  reasons  two  snubbo*  circuits  are 
included  across  the  power  diode.  Since  in  formulating  the  ACSL  simulations  it  is  assumed 
that  the  snubber  circuits  are  prqperly  designed  and  the  switching  behave  in  an  ideal  manner, 
these  circuit  elements  are  omitted.  Resistors  R29  through  R31  and  R35  through  R37  are 
bleeder  resistors.  These  resistors  are  7.5  k  ohms  and  are  needed  for  2  reasons:  to  equalize 
the  voltage  distribution  between  series  capacitors  and  to  discharge  the  cqjacitor  static 
charges  after  power  shut  down.  The  bleeder  resistors  are  ignored  in  the  ACSL  simnlattons 
The  total  bleeder  resistance  is  equivalent  to  11.25  k  ohms.  This  is  negligible  in  parallel 
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combination  with  the  load  resistance,  since  the  load  resistance  must  remain  below  114.75 
ohms  for  the  converter  to  remain  in  the  continuous  mode  of  operation.  In  order  to  obtain 
100  kW  across  the  output  with  an  output  voltage  of  750  volts  the  load  resistance  must  equal 
5.625  ohms.  A  resistance  load  smaller  than  5.625  ohms  will  require  the  converter  to  supply 

higher  than  rated  power. 

The  meter  devices  and  LED  indicators  at  the  input  and  the  output  terminals  of  Figure 
3.4  are  ignored  because  they  draw  negligible  amounts  of  current.  The  LRC  input  filter  is 
designed  fw  a  12-pulse  rectifier  configuration.  The  ACSL  simulations  including  active 
sources,  discussed  in  Chapter  V,  are  implemented  with  a  6-pulse  rectifier.  The  appropriate 
LRC  input  filter  design  is  documented  in  Chapter  V.  Figure  3.5  iUustrates  the  PSSCM 
configuration  based  on  the  given  design  and  simplifying  assumptions.  By  (3.19),  the  cutoff 
fiequency  for  the  converter  ouqiut  filter  is  found  to  be  85  Hz,  which  is  well  below  the 
intended  switching  frequency  of  5  kHz. 

f  =  _ ^ -  (3.19) 

^  2*PIVL1*C1 


Figure  3.5:  PSSCM  Configuration 
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IV.  VOLTAGE  CONTROLLERS  FOR  PSSCM 


Several  strategies  exist  to  control  the  duty  cycle  of  a  converter.  The  simplest  control 
implementation  is  a  fixed  duty  cycle.  The  fixed  duty  cycle  does  not,  however,  yield 
acceptable  dynamic  response  in  the  presence  of  varying  loads  or  varying  source  voltages. 
A  slightly  more  sophisticate  controller  is  the  Direct  Duty  Ratio  controller.  This  is  a  simple 
controller  in  which  the  output  voltage  is  compared  to  a  desired  voltage  and  multiplied  by  a 
gain  constant  to  obtain  the  duty  cycle.  Though  simple  to  implement  and  design,  the 
resultant  settling  time  and  overshoot  of  the  ou^ut  voltage  is  still  imacceptable.  This  chapter 
includes  a  discussion  covering  the  performance,  characteristics  and  function  of  three 
reliable  PSSCM  controllers:  the  PID,  the  Function  Control  and  the  standard  Multi-Loop 
controller. 

In  assessing  the  implementation  of  these  control  strategies,  attention  will  be  focused 
on  the  attractive  aspects  of  each.  Three  desirable  features  common  to  any  control  strategy 
are  given  as  follows  [4]: 

1.  The  output  voltage  of  the  regulator  remains  unchanged  even  though  there  are 
disturbances  from  either  the  supply  voltage  load  current. 

2.  The  closed-loop  equation  should  predict  the  performance  of  the  regulator. 

3.  The  control  circuit  should  be  simple  and  flexible. 

A.  THE  PID  VOLTAGE  CONTROLLER 

The  PID  voltage  controller  is  comprised  of  a  proportional,  integral  and  derivative  gain 
constant  and  hence  the  name  PID.  The  equation  governing  the  change  in  duty  cycle  about 
the  nominal  value  is  given  by  (4. 1).  It  is  a  function  of  the  change  in  ouq^ut  voltage  from  the 
desired  reference  voltage,  the  integral  of  this  change,  and  the  derivative  of  this  change 
times  various  gains.  There  are  no  feedforward  components  from  the  input  voltage  or 
feedback  components  from  the  inductor  or  ouq)ut  current.  If  the  capacitor  is  selected  large 
enough,  the  derivative  term  contribution  should  be  small. 
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(4.1) 


Ad=-K^p  (AVci)  -KpAVci  -  KiJAVcidt 
The  first  step  in  deriving  the  gain  constants  is  to  utilize  equations  (3.2)-(3.3)  to  derive 
the  state-space  averaged  model.  Equations  (3.2)  and  (3.3)  model  the  inductor  current 
dynamics  when  the  switch  is  in  the  closed  and  open  positions,  respectively.  Multiplying 
through  by  the  appropriate  time  duration  that  each  equation  is  vaUd,  d  or  d’,  equations  (4.2) 

and  (4.3)  are  derived. 

p(iLi)*d=  d*(El-Vci)/Ll  (4.2) 

p(iLi)*d’=  d’*(-Vci  /  LI)  (4.3) 

Since  d  4d’=l,  the  addition  of  (4.2)  and  (4.3)  yields  the  averaged  state-space  representa¬ 
tion  given  in  (4.4). 

pdu)  =(n)(<i»E-<fa) 

where  the  hat,  “  "  ",  denotes  averaged  quantities. 

The  output  voltage  is  made  avaUable  to  the  controUer  by  sensing  the  voltage  across  the 
cr^acitor.  The  equation  describing  the  dynamic  behavior  of  Vci  is  given  by 

Vci=(l/Cl)*Jic(dt)  (4-5) 

The  capacH^w'  current  is  equal  to  the  difference  between  the  inductor  current  and  the 
resistive  load  current.  The  current  through  the  load  is  equal  to  the  voltage  ao-oss  the  load 
divided  by  the  load  resistance. 

ic=iu-Vci/Rl 

The  averaged  voltage  state  equation  is  obtained  by  combining  (4.5)  and  (4.6)  and 
taking  the  derivative.  This  yields 


where  Vcj  is  the  average  c^acitor  voltage  and  Ili  is  the  average  inductor  current. 
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R1*C1  Cl 

li  0  ^ 

LI  LI 
aa  bb  cc 


(4.15) 


The  roots  of  A-SI  are  tamed  the  eigenvalues  (or  poles)  of  the  system  and  dictate  the 
stability  and  transient  behavior  about  an  operating  point.  Since  A  is  a  three-by-three  matrix, 
the  system  has  3  poles  in  the  transfer  function.  For  the  system  considered,  these  poles  are 


foimd  from 


—li _ S  —  0 

R1*C1  Cl 

A-IS=  _i  E 

LI  LI 

aa  bb  cc-S 


(4.16) 


The  characteristic  equation  of  this  third-order  system  is  obtained  from  the  determinant  of 
(4.16) 

S^  +  bS^  +  cS+d  (4.17) 


Where  the  coefficients  are 


^  R1*C1 


_  1  CC  El*l^ 

^  “  lI^  "  R1*C1  ”  LI 

El*aa  cc  E*bb 
^  ~  C1*L1  C1*L1  L1*R1*C1 


(4.18) 


(4.19) 


(4.20) 


Since  the  coefficients  of  the  characteristic  equation  are  a  function  of  the  gain 
constants,  the  selection  of  the  gain  values  will  affect  the  stability  of  the  system.  The  gain 
constants  are  obtained  by  simultaneously  solving(4.18)  through  (4.20) 


Kd  = 


L1*C1 

El 


)(^  Rl*Cl) 


(4.21) 
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■.-mi-ufel) 


(4.22) 

(4.23) 


The  required  controller  gains  can  thus  be  specified  once  a  desired  characteristic 


equation  is  selected  and  the  system  paramet^  are  known. 


B.  THE  FUNCTION  CONTROL  VOLTAGE  CONTROLLER 

The  information  used  to  implement  the  Function  Controller  is  found  in  [4].  The  output 
voltage,  Vc,  is  a  function  of  the  input  voltage,  E,  the  duty  cycle,  d,  and  other  intermediate 
variables  of  the  switch  converter  known  as  x: 

Vc=f(EAd)  (4.24) 

The  duty  ratio  is  a  function  of  the  ouq)ut  voltage,  reference  voltage,  Vref,  and  a 
combination  of  variables  of  the  switching  converter  known  as  y: 

d=f(y,Vc,Vref)  (4.25) 

The  key  element  of  the  derivation  of  the  duty  cycle  for  the  Function  Control  is  to  have 
the  output  voltage,Vc,  proportional  to  only  the  reference  voltage, Vref.  In  order  to 
accomplish  this,  the  duty  cycle  equations  for  the  operating  point  and  control  circuit  must 
be  equivalent.  The  operating  point  equation  and  the  control  circuit  equation  respectively 
can  be  obtained  from  (4.24): 

d<^=h(EAVc)  (4.26) 

d5e=h(E,x,K(Vref-Vc))  (4.27) 

Setting  (4.26)  and  (4.27)  equivalent  gives  Vc=K(Vref  -Vc),  which  in  turn  yields, 

Vc=  j^Vref  (4.28) 

where  for  large  K,  Vc  is  approximately  equal  to  Vref.  The  ouq)ut  voltage  is  now  a  function 
of  only  the  reference  voltage.  By  Kirchhoff ’s  voltage  law  the  following  equations  can  be 
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derived  from  the  low-frequency  averaged  equivalent  circuit  of  a  Buck  converta-  shown  in 
Figure  4.1.  In  this  representation  the  switch  is  modeled  by  the  average  current  assumed  to 
flow  through  it  while  the  diode  is  modeled  is  modeled  by  the  assumed  average  voltage 
across  it.  At  low  frequencies  the  inductor  voltage  and  the  capacitor  current  is  negligible. 


Figure  4.1:  Low-Frequency  Equivalent  Buck  Converter 


The  capacitcff  voltage  is  given  by 

Vc=d*E  -  L*p(iL)  (4-29) 

By  rearranging  (4.29)  the  duty  cycle  can  be  derived  for  the  operating  point: 

d=(Vc  +  L*p(iL))/E  (4.30) 

By  substituting  the  control  parameter  equation  (4.28)  into  equation  (4.30),  the 
duty  cycle  is  obtained  as 

d=(K(Vref-Vc)  +  L*p(iL))/E  (4.31) 

The  duty  cycle  for  the  Function  Control  is  a  function  of  the  difference  in  ouq^ut 
voltage  and  one  desired  reference  voltage,  the  derivative  of  the  current  through  the  inductor 
and  the  input  voltage.  In  particular  as  seen  in  (4.31),  the  duty  cycle  is  obtained  by  dividing 
through  by  E,  the  feedforward  input  voltage.  Division  circuits  are  usually  more 
complicated  to  implement  than  normal  gain  circuits  because  of  the  extra  hardware 
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involved.  An  analog  divider  followed  by  a  multiplier,  a  ratio-comparator,  or  a  digital 
interface  circuit  can  be  used. 


C.  THE  MULTI-LOOP  VOLTAGE  CONTROLLER 

The  multi-loop  controller  is  a  multi-variable  feedback  control  in  which  the  ou^ut 
voltage  and  the  inductor  current  are  both  used  to  adjust  the  duty  cycle.The  state  variable 
components  of  cuirent  and  voltage  are  both  incorporated  into  the  duty  cycle  perturbation. 
This  controller  has  no  voltage  derivative  component  nor  input  voltage  feedforward 
component.  The  variation  in  duty  cycle  is  given  by 

Ad=-hjAiLi  -  h^AVci  -  h^JAVcidt  (4.32) 

The  analysis  of  the  multi-loop  controller  is  similar  to  that  of  the  PID  voltage  controller. 
Upon  taking  the  doivative  of  (4.32)  and  substituting  (4.8)  and  (4.9),  the  following  state 
equation  is  derived 

+  (4.33) 

The  linearized  system  in  matrix  form  is  written  as 
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This  system  of  equations  is  once  again  in  the  standard  state  space  form,  px-B  x,  where 


X  is  the  state  vector  and  B  is  the  system  matrix 


—ll - L  0 

R1*C1  Cl 


0  4 

LI  LI 
am  bm  cni 

The  characteristic  equation  may  be  found  from 


(4.38) 


B-IS= 


— - S  —  0 

R1*C1  Cl 

-S  ^ 
LI  ^  LI 

am  bm  cm  -  S 


(4.39) 


As  discussed  before,  the  poles  of  any  third-order  system  are  determined  from  the 
characteristic  equation.  Using  standard  linear  algebra  techmques,  the  charactoistic 
equation  of  the  muli-loop  system  may  be  ejqjlicitly  written  as 


3  2 

S  +bmS  +cmS+dm 


where  the  coefficients  are  found  to  be 


(4.40) 


bm  = 


cm  = 


Ll+Rl*Cl*E*hi 
'  L1*C1*R1 

E*hj  +  Rl+R1*E*\ 
L1*C1*R1 


(4.41) 


(4.42) 


dm  = 


L1*C1 


(4.43) 


The  gain  constants  are  obtained  by  simultaneously  solving  (4.41)  through  (4.43).  This 


yields 


Ll*Cl*Rl(bm)  -  LI 
E1*R1*C1 


(4.44) 
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(4.45) 


.  ^  Ll*Cl(dm ) 

“  El 

Ll*Cl*Rl(cm)  -El(hi) 
__ 


(4.46) 


D.  CONTROLLER  COMPARISON  AND  IMPLEMENTATION  FOR  PSSCM 
The  PSSCM  specifications  determined  by  the  Power  System  Group  are  shown  in 
Table  4.1  [3].  The  PSSCM  configuration  is  based  on  the  assumptions  introduced  in  Ch^ter 
m  and  the  specifications  listed  in  Table  4.1.  The  simplified  circuit  is  shown  in  Figure  4.2. 

The  response  time  and  poformance  for  aU  three  controllers  discussed  in  this 
ch^ter  are  a  function  of  the  control  gains.  Based  on  a  desired  characteristic  polynomial, 
e}q)ressions  for  these  gains  were  derived  in  closed  form.  The  multi-loop  gnins  and  the  PID 
voltage  controllers  gains  determine  the  closed-loq)  pole  locations  of  the  system.  The 
characteristic  equations  of  these  two  systems  are  third  order.  There  are  several  techniques 
for  specifying  the  pole  locations.  The  following  pole  design  gave  the  best  results.  The 
placement  of  the  non-dominant  pole  should  be  approximately  one  decade  less  than  the 
switching  fi-equency.  The  switching  frequency  in  radians/sec  is  31.4  krad/sec,  therefore  the 
first  pole  win  be  fixed  at  3  krad/sec.  To  avoid  large  overshoots,  the  dominant  poles  are 
situated  in  the  over  damped  zone.  The  second  and  third  pole  are  complex  conjugates  and 
are  positioned  at  300  +/-  200j.  The  pole  placements  are  illustrated  in  Figure  4.3.  The  poles 
listed  in  Figure  4.3  lead  to  the  characteristic  equation 

S^  + 36008^ +  1.93e6S  +  3.9e8  (4.47) 
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Figure  4.2:  Simplified  PSSCM  Configuration 


TABLE  4.1:  PSSCM  PARAMETERS 


Parameters 

Values 

power 

lOOkw 

El 

850v+/- 

25v 

Vci 

750v 

R1 

5.625  ohms 

Cl 

2600  uf 

LI 

1.35  mH 

frequency 

5  khz 

nominal  D 

.8824 
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Figure  4.3:  PSSCM  Pole  Placement 


1.  Implementation  of  PH)  Controller  in  ACSL 

Using  the  coefficients  in  (4.47)  and  the  parameters  listed  in  Table  4.1,  (4.21) 
through  (4.23)  are  evaluated  and  provide  the  gains 


=  1.458e  -  5 

(4.48) 

Kp  =  .0068 

(4.49) 

Kj  =  1.6105 

(4.50) 

The  MATLAB  file  fw  calculating  these  gains  is  listed  in  ^pendix  A.  The  ACSL  com¬ 
mand  for  obtaining  the  duty  cycle  for  this  controller  is  also  provided  in  Appendix  B.  The 
derivative  component  of  the  control  is  implemented  by  a  differentiation  over  a  first  order 
lag. 

2.  Implemwitation  of  Control  Function  Controller  in  ACSL 

The  Function  Controller  was  designed  with  an  arbitrary  gain,  K=50.  Therefore 
using  (4.28),Vref  is  found  to  be  765  volts.  The  ACSL  command  fw  obtaining  the  duty 
cycle  for  with  this  controller  is  also  listed  in  Appendix  B. 

3.  Implementation  of  Multi-Loop  Controller  in  ACSL 

The  multi-loop  gains  can  be  obtained  by  stepping  through  the  same  procedures  as 
given  for  the  PID  controller.  Substituting  the  coefficients  of  (4.47)  and  the  parameters  listed 
in  Table  4.1  into  (4.44)  through  (4.46),  the  following  values  for  gains  are  obtained: 


hj  =  0.0056 

(4.51) 

h„  =  1.6105 

(4.52) 

h^  =  1.6105 

(4.53) 

The  MATLAB  file  for  calculating  these  gains  is  also  listed  in  Appendix  A  and  the  ACSL 
command  for  determining  the  duty  cycle  is  given  in  Appendix  B. 
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4.  Controller  Simulation  Comparisons 

To  contrast  the  performance  of  the  controllers  discussed  in  the  previous  sections, 
an  identical  ACSL  simulation  was  structured  and  executed  for  each.  In  the  program,  the 
current  is  obtained  by  integrating  the  derivative  portion  of  the  current  given  by  (3.3)  and 
(3.7),  respectively,  for  the  switch-on  and  switch-off  durations.  The  voltage  is  obtained  by 
integrating  the  current  in  the  Buck  capacitor  and  dividing  by  the  (xqjacitor  value  of  2600uf  . 
The  detailed  models  of  the  controllers  and  Buck  converter  are  used  in  the  simulations.  The 
simulation  runs  for  2.0  seconds  with  an  input  voltage  drop  at  1.0  second  and  a  resistance 
increase  at  13  seconds.  At  1.0  seconds  the  input  voltage  to  the  converter  is  reduced  from 
850  volts  to  800  volts.  This  voltage  drop  represents  a  change  of  approximately  6  percent. 
At  13  seconds  the  resistance  is  increased  by  a  factor  of  four  from  5.625  ohms  to  223  ohms. 
Figures  4.4  through  4.6  illustrate  the  results.  The  ou^ut  voltage  and  inductor  current  and 
the  duty  cycle  are  presented  in  each  figure  for  comparison  purposes.  The  start-up  transition 
prior  to  03  seconds  is  omitted,  so  that  the  resistor  load  changes  and  voltage  source  change 
can  be  enhanced  in  the  illustrations. 

The  simulation  results  show  the  Function  Control  controller  yields  the  least 
acceptable  transient  behavior.  Though  it  performs  fairly  well,  it  does  not  have  the  quick 
response  time  desired.  It  takes  sevo-al  tenths  of  a  second  to  stabilize  following  a  transient. 
The  PID  voltage  controller  and  the  muli-loop  controller  have  a  response  timft 
q>proximately  100  times  faster.  Lacking  an  integrator  in  its  structure,  the  Fimction 
controller  is  not  capable  of  guaranteeing  zero  steady-state  error  in  the  output  voltage 
response.  The  best  feature  of  the  Function  Control  controller  is  that  the  ouq>ut  voltage  peak 
spikes  are  limited  in  magnitude.  The  maximum,  a  spike  of  iq^proximately  6.5  volts,  occurs 
when  the  load  resistance  is  increased.  This  spike,  howevo",  is  only  a  0.87  p^cent  change 
from  the  desired  ou^ut  voltage  level. 

The  graphs  show  that  the  PID  voltage  controller  and  the  multi-lo<^  controller 
perform  remarkably  similar.  Close  analysis  reveals  that  the  PID  voltage  controller  is  in  fact 
similar  to  the  multi-loop  controller  because  the  capacitor  ciirrent  is  related  to  the 
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voltage  and  the  inductor  current  is  related  to  the  inductor  voltage  by  the  following 
equations: 

ici=Cl*p(Vci)  (4-54) 

‘u  =  (n)pu<it 

The  PID  voltage  controller  q)peaied  to  respond  to  the  source  voltage  drop  with 
Uttle  noticeable  effect  on  the  ouq)ut  voltage,  while  the  multi-loop  controller  had  a  sUght 
ripple  of  2  volts  in  the  output  voltage.  This  represents  a  variation  of  only  0.267  percent, 
which  for  all  practical  purposes  is  also  negligible.  The  load  increase  was  handled  slightly 
better  by  the  multi-loop  controller  than  the  PID  voltage  controller  because  the  PID  voltage 
contrnllftr  had  the  inductOT  current  spike  down  close  to  the  discontinuous  mode  of 
operation.  The  cuirent  spike  of  the  multi-loop  controller  is  approximately  12  amps,  which 
is  less  than  half  that  exp«ienced  with  the  PID  controller.  Both  controllers  had 
i^proximately  equal  ouq)ut  voltage  responses  with  a  maximum  peak  spike  of  8  volts  which 
is  only  a  one  percent  deviation.  It  can  be  seen  that  both  of  these  converters  compensated  for 
the  changes  in  load  and  input  voltage  within  milli-seconds  and  are  reliable  at  maintaining 
a  desired  ouq>ut  voltage.  Realistically,  since  load  changes  are  more  common  to  the 
proposed  system  than  step  voltage  source  disturbances,  the  performance  of  the  multi-loop 
/.ontmllftT  is  selected  as  the  candidate  fw  future  study  in  simulations  proposed  in  Ch^ters 
V  and  VI. 
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Figure  4.5:  Function  Control  Controller  Response 


Figure  4.6:  Multi-Loop  Controller  Response 


V.  BASIC  PSSCM  SIMULATIONS 


In  this  chapter,  simulation  results  are  presented  for  various  PSSCM  studies 
incorporating  the  multi-loop  control.  The  studies  include  step  changes  in  the  load 
resistance,  a  soft  source  voltage,  and  a  vector-controlled  induction  marhitift  start-up.  A  set 
of  representative  resrilts  is  presented  together  with  aritical  analysis. 

A  CONSTANT  SOURCEA^ARIABLE  LOAD  SIMULATIONS 

In  this  section,  the  main  focus  is  to  uncover  the  limits  and  1  ranges  of  current,  voltage, 
and  power  for  various  resistive  loads  on  a  PSSCM.  The  ideal  limits  and  ranges  are 
calculated  from  an  average- value  Buck  converter  model,  while  the  simnlfltpd  performance 
is  obtained  from  the  detailed  model  with  the  multi-loop  controller. 

1.  Ideal  Buck  Converter  Configuration 

This  section  documents  the  simulation  results  from  an  ACSL  program  that 
models  a  fixed  input  voltage  to  the  PSSCM  and  maintains  a  constant  ou^ut  voltage.  The 
ou^ut  voltage  is  set  at  750  and  the  input  voltage  is  set  at  850.  The  duty  cycle  varies 
depending  on  whidi  mode  of  opo'ation  the  converter  is  operating  in.  The  two  modes  of 
operation,  as  discussed  in  Chapter  HI,  are  continuous  current  and  discontinuous  current. 
The  load  resistance  value  will  dictate  the  operating  mode  since  the  frequency  and  the 
inductor  parameters  are  already  fixed.  Equation  (3.13)  establishes  the  boundary  condition 
for  the  two  modes.  Rearranging,  an  expression  for  the  critical  resistance  is  given  by 

R.  =  (5,„ 

Using  (5.1),  the  oitical  value  of  resistance  is  found  to  be  114.75  ohms.  Any 
resistive  load  less  than  that  will  keep  the  converter  in  continuous  mode,  while  a  resistive 
load  above  that  level  will  cause  the  converter  to  operate  in  discontinuous  mode. 
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Figure  5.1  shows  the  inductor  current  plotted  for  4.0  ms.  In  this  study  the 
resistance  value  is  changed  from  5.625  ohms  to  114.75  ohms  at  1.2  ms  and  then  to  150 
ohms  at  2.4  ms.  As  shown  in  Figure  5.1,  the  inductor  current  goes  from  continuous 
operation  with  a  minimum  current  of  126.8  amps  to  discontinuous  mode  with  zero 
minimum  current.  While  operating  with  the  critical  value  of  load  resistance,  the  current 
only  instantaneously  reaches  zero,  while  after  2.4  ms,  operating  with  150  ohm  resistance, 
the  current  is  zero  fw  a  duration  of  20.6  us.  This  is  approximately  10  percent  of  the  duty 

cycle. 

Figure  5.2  shows  the  inductor  current  and  load  power  for  a  series  of  simulated 
step  changes  in  load  resistance.  The  load  resistance  varies  from  5.625  to  150  ohms.  The 
step  changes  include  a  doubling  of  resistance  from  5.625  to  11.25,  then  to  25  ohms,  then 
increases  in  steps  of  25  ohms  to  150  ohms.  As  the  load  resistance  doubles,  the  current  and 
power  are  reduced  in  half.  Table  5.1  contains  the  ACSL  computed  results  of  the  minimum 
and  maximum  current,  the  load  power  and  the  mode  of  operation.  The  table  also  contains 
the  duty  cycle,  D,  and  an  additional  parameter  D2.  D2  indicates  the  point  in  a  given  cycle 
at  which  the  current  goes  discontinuous. 
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Figure  5.1:  Ili  Continuous  and  Discontinuous  Mode 
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Note  that  in  discontinuous  mode  the  duty  cycle  is  different  then  in  continuous 
operation.  When  the  switch  is  closed  during  continuous  operation,  the  current  is  linearly 
reduced,  but  in  discontinuous  operation  the  current  goes  non-linear  at  D2.  The  ACSL 
program  for  the  ideal  Buck  with  constant  input  voltage  and  steady-state  output  voltage  is 
included  in  ^pendix  C. 
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Figure  5.2:  Step  Changes  in  Resistance-Constant  InpuVOu^ut  Voltage 


TABLE  5.1:  Data  for  Different  Load  Resistance  Conditions 


Rl(ohms) 

^mtn 

(amps) 

tmax 

(amps) 

Power 

(Watts) 

operation 

mode 

D 

D2 

5.625 

126.80 

139.87 

100k 

cm 

.8823 

NA 

11.25 

60.13 

73.20 

50k 

cm 

.8823 

NA 

25 

23.46 

36.54 

22.5k 

cm 

.8823 

NA 

50 

8.46 

21.54 

11.25k 

cm 

.8823 

NA 

75 

3.46 

16.54 

7.50k 

cm 

.8823 

NA 

100 

0.96 

r  14.04 

^  5.63k 

cm 

.8823 

^NA 

125 

0.0 

12.53 

4.50k 

dcm 

.845 

.958 

150 

0.0 

11.43 

3.75k 

don 

.772 

.875 

2.  Basic  Buck  Converter  with  Related  Output  Voltage 

In  this  section,  the  simulation  results  of  the  Buck  converts  with  stepped  load 
resistance  are  presented.  The  input  voltage  source  to  the  converter  is  fixed  at  850  volts.  The 
multi-loop  controUer  introduced  in  Chapter  IV  is  used  to  establish  the  duty  cycle  of  the 
converter.  The  stabiUty  of  the  output  voltage  and  the  behavior  of  the  inductor  current  are 
the  main  focus  of  these  studies.  The  ACSL  program  for  the  simulations  in  this  g^ytion  is 
included  in  Appendix  D.  As  ^^ical  PSSCM  loads  will  guarantee  continuous  current 
operation,  the  stepped  resistance  values  will  be  kept  below  the  critical  value.  The  ou^ut 
voltage  is  determined  by  integrating  over  the  c^acitor  current  and  dividing  by  the 
capacitance 

(5.2) 

where, 

ici  =  iLl-iRi  (5.3) 
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Figure  5.3  illustrates  the  simulation  results  for  the  ou^ut  voltage  and  inductor 
C^t  for  stop  ohangos  in  load  from  5.625  to  100  ohms.  The  output  voltage  is  weU 
regulated  with  only  smaU  voltage  transients.  The  largest  transient  occurs  when  the  load  m 
increased  from  5.625  ohms  to  11.25  duns  with  a  voltage  spiie  of  approximately  5  volts. 
Hus  transient  is  deemed  negligible,  representing  only  a  0.66  percent  vanatton  m  the 
nominal  750  volt  output  voltage.  As  the  resistance  approximately  doubles  from  1 1.2  to  25, 
25  to  50.  and  50  to  100  ohms,  the  magnitude  of  the  output  voltage  transients  decrcase  by 

approximately  50  percent  during  each  step  change. 

Figure  5.4  illustrates  step  changes  in  load  resistance  from  5.625  to  25, 25  to  100, 
100  to25.  and  25  to  5.625  ohms.  Each  step  change  is  four  times  the  previous  value  for 
iucreasing  increments  and  one-fourth  for  decreasing  increments.  A  small  voltage  spike 
occurs  at  0.2  seconds  when  the  resistance  changes  from  5.625  to  25  ohms.  The  spike  at  this 
point  is  approximately  8  volts  which  is  one  percent  of  the  nominal  output  voltage.  The 
laigest  spike  occurs  at  0.8  seconds  when  the  load  resistance  decreases  from  25  to  5.625 
ohms.  The  transient  spike  is  approximately  25  volts,  which  represents  a  3.3  percent 
variation  from  the  nominal  output  voltage.  At  this  point  there  is  also  a  noticeable  35  amp 

current  spike  in  the  inductor  current  This  overshoot  is  approximately  23  percent 

Figure  5.5  illustrates  larger  increases  in  load  step  changes.  The  changes  are  from 
5.625  to  50. 50  to  100, 100  to  50,  and  50  hack  to  5.625  ohms.  There  is  an  approximate  10 
volt  spike  «  the  load  resistance  increases  from  5.625  to  50  ohms.  This  is  roughly  a  1.3 
percent  variation  from  the  output  voltage  which  is  only  sUghtly  larger  than  the  8  volt  spike 
from  5 .625  to  25  ohms.  However,  there  is  enough  of  a  decrease  in  current  to  cause  the 
inductor  current  to  instantaneously  go  discontinuous.  Again,  the  largest  spike  occurs  at  0.8 
seconds  as  the  load  resistance  decreases  from  50  to  5.625  ohms.  The  spike  is  about  30.6 
volts  which  is  4  percent  of  the  nominal  output  voltage.  At  this  point  the  current  overshoot 

is  approximately  45  amps  or  35  percent. 
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Figure  5.6  illustrates  step  changes  in  load  resistance  from  5.625  to  100  ohms  and 
then  back  down  to  5.625  ohms.  As  the  resistance  is  ino’eased,  the  inductor  current  goes 
discontinuous  fm  a  few  milli-seconds  and  a  voltage  spike  of  approximately  12.5  volts  is 
obsoved.  At  0.6  seconds  when  the  resistance  drops  back  down,  the  current  overshoot  is  the 
greatest,  q>proximately  35  amps  or  38  percent  The  largest  voltage  spike,  35  volts,  also 
occurs  at  this  point  Though  representing  the  largest  spike,  it  is  still  only  5  percent  of  the 
nominal  ou^ut  voltage.  The  muti-loop  controller  adequately  regulates  the  ouqjut  voltage 
for  large  and  small  changes  in  load  resistance. 
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Figure  5.4:  Resistive  Load  Step  Change  for  Constant  Source  (Initial  step  25  ohms) 
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Figure  5.5;  Resistive  Load  Step  Change  for  Constant  Source  (Initial  step  50  ohms) 
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Figure  5.6:  Constant  Source  Full  Load  Step  Change 


B.  ACTIVE  SOURCEA^ARIABLE  LOAD  SIMULATIONS 

In  this  section,  the  operation  and  implementation  of  an  active  source  is  introduced. 
This  source  is  then  connected  to  the  PSSCM  and  operation  under  various  load  conditions 
is  investigated.  In  particular,  studies  are  presented  for  changes  in  load  resistance  and  the 
{^plication  of  a  inverter^duction  motor  load. 

1.  Design  of  Active  Source 

The  active  source  configuration  has  a  3-phase  voltage  source,  a  rectifier  stage  and 
a  low-pass  filter  stage  as  depicted  in  Figure  5.7.  This  is  an  “active”  source  because  the 
outi)ut  voltage,  E,  will  vary  as  the  current  drawn  from  the  rectifier  varies. 


The  3-phase  line-to-neutral  voltage  source  can  be  converted  to  a  3-phase  line-to- 
linft  voltage  source  by  scaling  the  amplitude  by  the  square  root  of  three  and  shifting  the 
phases  by  thirty  degrees.  The  line-to-neutral  voltages  are  denoted  by  Vag,  Vbg  and  Vcg 
while  the  line-to-line  voltages  are  given  by  Vab,  Vbc  and  Vca.  The  3  -phase  bridge  rectifier 
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convols  the  3-phase  AC  voltages  to  a  DC  voltage  with  hannonics  at  multiples  of  six  times 
the  source  firequency.  The  average  voltage  out  is  0.955  Vmax(abc).  The  rectifier  is  referred 
to  as  a  six-pulse  converter  because,  for  continuous  i^tho'e  are  6  segments  to  eadi  period  in 
the  ouqjut  waveform.  A  detailed  discussion  of  the  operation  of  the  rectifier  may  be  found 
in  Chapter  5  of  Fisher  Power  Electronics  [5].  The  RLC  low-pass  fUto'  reduces  the  ripple 
voltage  to  a  smooth  usable  DC  voltage.  The  fundamental  frequency  of  the  ripple  is  360  Hz 
far  a  3-phase  6-pulse  bridge  rectifier.  For  a  12-pulse  rectifier  the  ripple  frequency  is  720 
Hz.  Equation  (3.18)  dictates  the  necessary  constraints  required  to  specify  a  suitable  cutoff 
frequency  to  eliminate  the  harmonics  caused  by  the  rectifio:.  The  filte-  design  used  in  the 
active  source  has  a  10.3  Hz  cut-off  frequency  which  is  well  below  the  harmonics  caused  by 
the  rectifier.  The  design  values  for  the  filter  impedances  X^,  X^and  Rfaic  4.725, 0.1332 

and  0.236  ohms  respectively,  ^pendix  E  contains  a  listing  of  the  ACSL  program  for  the 
active  source. 


2.  Variable  Resistive  Load  Simulations 

Simulation  studies  identical  to  those  presented  for  the  fixed  source  in  Section  3.B 
are  shown  in  Figures  5.7  through  5.10  using  the  active  source.  The  transient  overshoot 
magnitudes  in  the  ouqjut  voltage  and  inductor  current  are  similar  to  those  encoimtered  with 
the  constant  source.  This  holds  for  small  and  large  step  changes  in  load  resistance,  up  or 
down.  The  primary  difference  between  the  active  source  simulations  and  the  fixed  nnngfflnt 
source  simulations  is  that  the  active  source  takes  longer  to  adjust  to  the  changes  in  load.  In 
particular,  the  active  source  requires  a  longer  time  to  respond  as  the  load  resistance  is 
stepped  upward.  The  second  difference  is  that  the  active  source  has  a  slight  r4>ple  in  the 
output  voltage.  As  the  load  resistance  increases,  the  ripple  deca-eases.  At  the  minimum  load 
resistance  level,  the  peak-to-peak  ripple  is  only  about  2.5  volts.  This  is  a  0.33  percent 
variation  in  the  ouq>ut  voltage.  To  try  to  ino-ease  the  PSSCM  performance,  the  multi-loop 
controller  gain  constants  discussed  in  Chapter  IV  were  set  to  a  variable  mode;  that  is,  the 
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parameters  h^.  b^rndK  were  changed  on-the-fly  as  the  load  resistance  was  changed.  The 
results  of  these  simulations  showed  no  improvement  in  response  and  in  general, 
comparable  levels  of  performance  as  illustrated  in  the  simulations  with  constant  gam 

parameters. 
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Figure  5.10:  Analysis  of  Load  Step  from  5.625  to50  ohms  with  Active  Soiiix:e 
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3.  Induction  Motor  Simulations 

A  460  Volt  50  hp  induction  machine  model  was  incorporated  as  a  load  on  the 
Buck  converter  and  various  simulation  studies  were  conducted.  The  machine  parameters 
for  the  induction  motor  are  listed  in  Chapter  4  of  Analysis  of  Electric  Machinery  [6].  The 
rated  speed  is  1705  ipm  which  is  equal  to  178.55  rad  /sec.  The  460  volt  ac  power  required 
to  run  die  motor  must  be  generated  by  a  Ship  Service  Inverter  Module  (SSIM).  The  power 
to  the  SSIM  is  the  750  volt  dc  power  provide  by  the  PSSCM.  The  inverter  peak  phase 
voltage  is  2/it  of  the  750  volt  input  to  the  SSIM,  which  is  478  volts.  The  initial  SSIM  was 
configured  as  a  six-step  polyphase  inverter  system  [5].  Without  any  current  regulation, 
upon  start-up  this  system  drew  unacceptable  levels  of  current  from  the  PSSCM.  As  a  result, 
it  also  caused  the  inductor  current  to  fluxuate  between  the  continuous  and  discontinuous 
modes  of  operation. 

The  SSIM  configuration  employed  in  the  ACSL  simulfltinng  uses  a  hysteresis 
current  controller.  Chapter  2  of  Current  Control  ofVSI-PWM  Inverters  [7]  and  chapter  17 
of  Mohan  Power  Electronics  [8]  contain  discussions  on  the  qieration  of  the  hysteresis 
controller. 

The  hysteresis  controlled  PWM  scheme  is  inherently  current  limiting,  so  large 
inrush  currents  caused  by  the  induction  machine  are  controllable.  The  hysteresis  controller 
in  conjunction  with  a  vector  controller  changes  both  the  commanded  currents  and  the 
cotxunanded  current  frequency.  The  vector  controller,  by  fixing  all  of  the  rotor  flux  to  be 
along  the  d-axis  of  the  synchronous  ref^ence  frame,  provides  near  instantaneous  torque 
control.  The  commanded  current,  is  determined  by  the  vector  controller.  The  phase 
current,  ijs,  is  subtracted  from  the  command  current  in  order  to  obtain  the  current  error.  If 
the  current  oror  is  positive  than  the  current  is  reduced  by  closing  the  top  switch  of  the 
inverter  leg.  If  the  current  error  is  negative  than  the  current  is  increased  by  closing  the  other 
switch  in  the  inverter  leg.  These  switches  are  thus  toggled  on  and  off  in  what  is  termed  a 
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pulse-width-modulation  pattern.  The  same  principle  appUes  to  the  other  two  legs  of  the 
inverter. 

The  current  and  voltage  equations  used  in  the  ACSL  program  to  convert  between 
the  stationary  and  synchronous  reference  frame  are  derived  in  Krause,  Electromechanical 
Motion  Devices  [9].  Appendix  E  contains  the  ACSL  code  for  the  simulation  of  the  inverter 
and  hysteresis  controUer  in  conjunction  with  the  active  source.  Appendix  F  contains  the 
induf  tinn  machine  macro  which  is  required  by  the  ACSL  program  in  ^pendix  E. 

Figures  5.12  and  5.13  portray  the  simulation  results  of  the  PSSCM  operating  with 
an  induction  motor  and  variable  load  resistance.  In  this  simulation  the  load  resistor  is 
inweased  from  10  ohms  to  20  ohms  at  1.0  seconds,  from  20  to  50  ohms  at  2.0  seconds,  and 
50  fthms  to  100  ohms  at  3.0  seconds.  The  induction  motor  is  in  parallel  with  the  load 
resistance,  and  is  turned  on  at  0.05  seconds  and  reaches  full  speed  in  approximately  1.0 
second.  The  change  in  resistance  has  negligible  effect  on  the  ouqjut  voltage  which  would 
be  expected  because  the  parallel  combination  of  motor  resistance  with  the  Buck  resistor 
reduces  the  total  load  resistance.  Though  the  average  output  voltage  remains  at  a  steady  750 
volts  the  induction  motor  increases  the  ripple  effect  by  about  +/-  8  volts.  The  induction 
motor  appears  as  i^jproximately  a  400  ohms  average  resistance  at  turn  on  and  reduces  to 
approximately  30  ohms  at  the  1705  rpm  rated  speed.  Once  the  induction  motor  reaches 
rated  speed,  changp-s  in  the  load  resistance  spears  to  have  no  effect  on  the  motor  speed. 
Even  when  the  load  resistor  is  removed,  the  maximum  total  average  resistance  not  exceed 
30  ohms  at  rated  speed.  The  power  required  by  the  induction  motor  at  rated  speed  is  18,750 
watts.  Table  5.2  contains  the  average  value  results  of  the  simulations  in  Figures  5.12  and 

5.13. 

Figures  5.14  and  5.15  contain  the  simulation  results  of  an  induction  motor  in 
parallel  with  a  load  resistance  of  100  ohms.  The  motor  starts  at  0.05  seconds  and  the  study 
extends  to  4.0  seconds.  Before  the  induction  motor  is  put  on-line,  the  external  resistance  is 
equal  to  the  converter  load  resistor.  Under  these  conditions,  the  current  operates  near  the 
boundary  of  the  discontinuous  mode.  When  the  induction  motor  is  put-on  line  the  total 
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average  resistance  falls  to  about  80  ohms.  This  represents  the  parallel  combination  of  the 
average  400  ohm  resistance  introduced  by  the  induction  motor  and  the  100  ohm  load 
resistor.  Once  the  motcH'  reaches  rated  speed  the  average  parameter  values  are  the  same  as 
those  documented  in  Table  5.2  with  the  100  ohm  resistor. 


TABLE  5.2:  Effects  of  Resistance  Changes  on  Induction  Motor 


R 

WRM 

P. 

^mv 

^total 

Vc 

II 

ohms 

ohms 

ohms 

r.p.m. 

kw 

kw 

volts 

amps 

20 

30 

12 

1705 

18.75 

47.0 

750 

62.5 

50 

30 

19 

1705 

18.75 

30.0 

750 

40.0 

100 

30 

23 

1705 

18.75 

24.5 

750 

32.6 

55 


Figure  5.12:  Variable  Resistance  Simulation  with  Induction  Motor 


Figure  5.13:  Variable  Resistance  Simulation  with  Induction  Motor 
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Fiffure  5.14:  Induction  Motor  Simulation  with  R=100  ohms 
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Figure  5.15:  Induction  Motor  Simulation  with  R=100  nhms 
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VI.  ADVANCED  PSSCM  SIMULATIONS 


k  this  chfq>ter,  the  perfcHtnance,  characteristics  and  functionality  of  three  specially 
configured  dual  converter  systems  are  discussed  and  various  simulation  study  results  are 
presented. 

A.  DUAL  CONVERTER  PARALLEL  CONFIGURATION 

The  dual  converter  parallel  configuration  represents  two  separate  PSSCMs  regulating 
two  independent  zones  off  of  the  same  bus  as  illustrated  in  Figure  6.1.  The  bus  voltage  is 
supplied  by  the  active  source  discussed  in  Chapter  V.  The  topological  configurations  for 
each  of  the  PSSCMs  are  identical. 


L2 

converter  2 

Figure  6.1:  Dual  Converter  Parallel  Configuration 


61 


Figures  6.2  and  6.3  illustrate  simulation  study  results  for  both  converters  running  460 
volt  induction  motors.  Initially  both  converted  have  100  ohm  load  resistors.  At  1.0  second 
the  induction  motor  in  zone  1  is  put  on-line,  and  at  3.0  seconds  the  induction  motor  in  zone 
2  is  put  on-line.  At  the  time  the  induction  motOT  is  put  on  line,  the  associated  hystoesis 
inverter  is  also  activated.  The  ou^ut  voltages  for  both  converters  are  found  to  be  regulated 
at  750  volts  with  a  +/- 12  volt  variation.  Each  PSSCM  is  regulated  by  an  independent  multi¬ 
loop  controller.  The  average  inductor  currents  of  both  converters  are  near  discontinuous 
until  the  motor  is  put  on-line,  then  the  currents  increase  to  an  average  of  approximately  65 
amps  at  rated  speed.  Due  to  the  designed  ramp  up  the  of  commanded  rotor  speed,  there  is 
no  noticeable  effect  on  converterl  when  the  second  induction  motor  in  zone  2  is  put  on- 
line.There  is,  however,  a  noticeable  drop  in  the  supply  voltage.  The  bus  voltage  drops  by 
about  13  volts  when  each  induction  motor  is  started.  To  simulate  more  than  2  zones  a  more 
stable  active  source  should  be  used. 

Figures  6.4  and  6.5  iUustrate  the  inter-zonal  effects  of  zone  2  transients  leading  to 
disturbances  in  zone  1.  In  particular,  this  simulation  illustrates  an  induction  motor 
operating  in  steady  state  in  zone  1  while  a  large  change  in  load  resistance  occurs  in  zone  2. 
At  25  seconds  the  zone  2  load  resistance  is  stepped  up  from  5.625  ohms  to  100  ohms  and 
at  3.2  seconds  the  resistance  is  reduced  back  to  5.625  ohms.  When  the  resistance  is  stepped 
up  in  zone  2,  the  voltage  and  inductw  current  in  zone  1  has  no  noticeable  transients,  but 
when  the  resistance  is  stepped  back  down  to  5.625  ohms  there  is  some  noticeable  effects. 
The  ou^ut  voltage  in  zone  1  has  transients  of  +/-  40  volts  caused  by  the  large  power 
requirements  of  zone  2.  In  fact,  the  power  requirement  ina-eases  from  5.6  kw  to  100  kw  a 
factor  of  nearly  18  times.  The  power  supply  to  the  converters  experiences  large  transients 
of  approximately  100  volts  caused  by  both  the  increase  and  decrease  in  load  resistance  in 
zone  1. 

To  simulate  additional  zones  a  steam  turbine  with  a  three-phase  exciter  can  be  used. 
The  block  of  tiiis  active  source  is  illustrated  in  Figure  6.2.  ^pendix  G  contains 

the  ACSL  command  for  this  active  source  including  the  macros  required  for  simulation. 
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Single  zone  simulations  poframed  with  tins  active  source  were  found  to  be  equivalent  to 
those  presented  for  a  constant  fixed  source.  This  follows  because  of  the  large  capacity  of 
the  associated  turbine/generator  combination. 


Figure  6.2:  Steam  Turbine  with  3-Phase  Synchronous  Machine 
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Figure  6.6:  Inter-Zonal  Effects 


B.  DUAL  CONVERTER  SERIES  CONFIGURATION 

The  dual  converter  series  configuration  simulates  the  PSSCM  providing  a  regulated 
source  of 750  volts  to  a  second  SSCM  in  the  same  zone.  The  secondary  SSCM  converts  the 
750  volt  DC  source  to  500  volts  DC.  The  500  volt  conversion  was  arbitrarily  selected  for 
simulation  purposes.  This  section  presents  an  analysis  of  the  stability  and  performance  of 
the  primary  system  with  changes  in  load  resistance  on  the  secondary  system. 


switch  1  LI 


PSSCM  SSCM 


Figure  6.7:  Dual  Converter  Series  Configuration 

Figures  6.8  and  6.9  illustrate  simulations  of  the  primary  and  secondary  systems  as  the 
secondary  converter  is  put  on-line  and  the  load  resistance  is  varied.  The  input  voltage  to  the 
primary  converter  is  held  constant  at  850  volts  and  the  resistor,  Rl,  is  100  ohms.  Both 
converters  have  the  same  inductor  and  cf^acitor  component  values,  switching  frequency, 
and  multi-loop  controller  design.  At  0.2  seconds  the  secondary  converter  is  activated  with 
the  ouq)ut  voltage  starting  from  an  initial  condition  of  zero.  Because  of  the  large  transient 
power  requirements,  the  output  voltage  of  the  primary  converter  experiences  a  sudden  drop 
from  750  to  600  volts.  This  is  a  transient  of  150  volts,  which  corresponds  to  20.0  percent. 
The  inductor  current  of  the  primary  converter  spikes  up  with  an  overshoot  of  195  amps.  The 
power  required  is  50  kW  for  a  resistor  of  5  ohms  at  500  volts,  which  is  50  percent  of  the 
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maximum  power  specifications  for  a  SSCM.  TTie  primary  converter  powa*  requirement  is 
ino-eased  from  7.23  kW  for  a  100  ohm  resistor  to  68  kW. 

These  transients  are  not  tolerable  and  must  be  reduced  by  having  a  larger  secondary 
converter  resistance,  adding  a  dampening  network  and/or  ramping  the  tum-on  phase  of  the 
secondary  converter.  At  0.4  seconds  the  resistance  of  the  secondary  convoter  is  stepped  up 
from  5  ohms  to  40  ohms  and  at  0.7  seconds  the  resistor  is  deceased  back  down  to  5  ohms. 
A  value  of  40  ohms  was  used  as  the  maximum  resistance  because  at  41  ohms  the  current 
goes  discontinuous.  The  changes  in  load  resistance  of  the  secondary  converter  have  little 
effect  on  the  primary  converter  ouq>ut  voltage.  Figures  6.8  and  6.9  illustrate  the  duty  cycles 
for  both  primary  and  secondary  controllers,  represented  by  D1  and  D2  respectively. 
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Figure  6.8;  Simulation  Effects  Caused  by  Secondary  SSCM 


C.  DUAL  CONVERTERS-INTERLACE  CONFIGURATION 

The  configuration  represents  an  attempt  to  provide  a  zone  with  a  more  stable 

voltage  source  by  integrating  two  partial  Buck  converter  topologies  together  in  a  PSSCM. 
The  design  uses  two  multi-loop  controllers  to  provide  the  voltage  regulation  for  a  given 
load  as  shown  in  Figure  6.10. 


controller  1  LI 


Figure  6.10:  Dual  Converter  Interlace  Configuration 

Ideally,  with  2  controllers  providing  regulation,  the  PSSCM  performance  should  be 
enhanced,  and  the  average  currents  through  the  switches  and  inductors  should  be  reduced. 
To  better  analyze  this  system  a  constant  source  of  850  volts  is  supplied  to  the  converters. 
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Several  simulations  were  performed  in  an  attempt  to  enhance  the  performance  of  the 
PSSCM.  Depending  on  the  simulation  performed,  at  least  one  of  the  following  parameters 
was  changed:  the  switching  frequency,  the  inductor  values,  the  multi-loop  contoller 
parameters  and/or  the  nominal  duty  cycle.  The  attempt  to  improve  the  performance  of  the 
PSSCM  was  unsuccessful.  In  most  cases  the  performance  was  degraded  iising  the 
interlaced  controllers. 
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m  CONCLUSION 

A.  OVERVIEW  OF  SHIPBOARD  MORE  -ELECTRIC 

The  design  of  the  Primary  Ship  Service  Converter  Module  (PSSCM)  controls  is  the 
key  to  providing  each  zone  with  regulated  and  stable  power.  A  representation  of  a  PSSCM, 
upon  employing  various  assumptions  and  reductions,  was  simulated  using  the  Advanced 
Continuous  Simulation  Language  (ACSL).  After  applying  several  assumptions,  the 
PSSCM  was  reduced  to  a  basic  DC-DC  Buck  converte*.  Fixed  constant  sources  and  active 
sources  were  specified  to  investigate  and  critique  the  performance  of  the  PSSCM. 

To  enhance  die  stability  of  the  PSSCM,  several  controllers  were  implemented  to 
control  the  duty  cycle.  From  the  three  voltage  controllers  analyzed,  the  multi-loop 
controller  yielded  the  most  desirable  performance.  The  multi-loop  controller  regulated  the 
ouq>ut  voltage  the  most  effectively  by  reducing  the  oscillations,  voltage  spikes  and 
response  times  during  transient^, 


B.  SUMMARY  OF  RESULTS 

The  basic  concepts  of  the  Buck  converter  were  introduced  and  applied  to  both  ainniiflte 
the  PSSCM  and  integrate  it  together  with  some  typical  zone  loads.  Several  studies  were 
analyzed  to  determine  the  stability  and  perfcxmance  of  the  PSSCM.  For  the  multi-loop 
controller,  it  has  been  shown  that  this  regulation  holds  sctoss  the  spectrum  of  large  and 
small  increases  and  decreases  in  load  resistance.  The  PSSCM  behavicH*  is  stable  for  both 
active  and  fixed  sources.  The  power  requirements  of  the  PSSCM  remain  within  toVrance 
as  long  as  the  total  output  converter  resistance  is  maintained  above  5.625  ohms. 
Appropriate  regulation  of  the  currents  required  by  active  loads,  such  as  induction  machines, 
must  be  maintained  fw  the  PSSCM  to  operate  properly.  It  was  shown  that  a  vector 
controlled  inverter  system  possessed  the  current  controlling  capabilities  required  for 
successful  operation  with  a  PSSCM. 
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The  «;»"iiiintinns  performed  to  study  the  effects  of  two  PSSCM  running  simultaneously 
from  the  same  power  source  illustrated  that  inter-zonal  effects  are  present  for  large  changes 
in  power  requirements.  The  effects  of  putting  secondary  converters  on-line  are  significant 
and  must  be  considered  to  avoid  intolerable  disturbances.  Though  for  less  demanding 
conditions,  it  has  been  shown  that  this  regulation  holds  across  the  spectrum  of  typical  load 
resistance  changes  and  induction  machine  operations.  To  increase  the  stabiUty  of  the 
converters  a  damping  resistor  or  a  RC  damping  network  cah  be  included.  A  damping 
network  would  enhance  the  transient  time  and  overshoot  performance.  Though  additional 
component  cost,  weight  and  space  lirmtation  should  be  considered. 

The  effort  to  increase  the  performance,  reliability  and  efficiency  of  the  Buck  converter 
by  interlacing  the  control  signals  implied  to  paralleled  convoters  was  unsuccessful. 


C.  FUTURE  RESEARCH  CONSIDERATIONS 

The  steam  turbine  active  source  discussed  in  Ch^ter  VI  can  be  implemented  together 
with  several  representative  zones.  This  larger  scale  analysis  would  give  a  better 
representation  and  further  insight  into  the  performance  and  stability  issues,  especially  inter¬ 
zonal  effects.  A  12-pulse  rectifier  system  can  be  designed  and  implemented  to  replace  the 
6-pulse  rectifier  system  discussed  in  Chapter  V.  The  non-ideal  effects  caused  by  the 
inductor  resistance  and  switch  resistance  can  be  incorpOTated  into  the  analysis  to  get  a  more 
comprehensive  model  of  the  SSCM. 

By  inrb^ding  in  the  analysis  and  simulations  snubber  circuits,  AC  filtering  c^acitors, 
bleeder  resistors  and  other  non-ideal  effects  discussed  in  Chq)ter  HI,  the  efficiency  of  the 
PSSCM  can  be  better  assessed. 
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APPENDIX  A.  MATLAB  PROGRAM  FOR  FINDING  VOLTAGE  CONTROLLER 

GAINS 


%  Determine  control  parameter  gains  for  controllers  for  PSSCM 
% 

f=5e3; 

T=l/f; 

E=850; 

Vc=750; 

D=Vc/E; 

d=0.0; 

dcap  =  (D  +  d)*T; 

Z=D*T; 

rl=5.625; 

r2=11.25; 

r=5.625; 

C— 2600e"'6j 

L=1.35e-3; 
w=2*pi*f; 
ml=-300  +j*200; 
m2=-300  -j*200; 
m3=-3000; 
m=[ml,m2,m3]; 

chm=poly(m);  %  results  gives:  1  3600  1.93e6  3.9e8 
% 

%  PID  voltage  controller  gains: 

kd=  (L*c/E)*(chm(2)  -  (l/(r*c)))  %  1.458e-5 

kp=  (L*c/E)*(chm(3)  -  l/(L*c))  %  .0068 
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ki=  (c*L/E)*(chm(4))  %  1.6105 


aa=  kp/(r*c)  -  ki  +  (-(l/((r*c)A2))  +  l/(L*c))*kd; 

bb=kd/(r*cA2)-kp/c; 

cc=  (-W(L*c))*kd; 

A=[-l/(r*c) ,  1/c  ,0;  -1/L ,  0 ,  E/L;  aa  ,bb,  cc]; 


%  multi-loop  controller  gains; 
hl=  (chm(2)*(L*r*c)  -  L)/(E*  r  *c)  %.0056 
hN=  chm(4)*(L*c)/E  %1.6105 

hV=  (chm(3)*(L*r*c)  -E*hl  -  r)/(E*r)%.(X)58 
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APPENDIX  B.  IMPLEMENTATION  OF  CONTROLLERS  IN  ACSL 


1. )  PH)  voltage  controller 
Kd=1.458e-5 
Ki=1.6105 
Kp=0.0068 

delVc=  (Vc-VcD) 

Vex  =  INTEG(delVc.0.0) 
fvh=INTEG(Vcy,0.0) 

Vcy  =  (Kd*delVc  -  INTEG(Vcy,0.0))/(TD) 
delD=  -Kp*delVc  -  Vcy  -Ki*Vcx 
D  =  BOUND(1.0e-5,1.0,  VcD/E  +  DelD) 

2. )  Function  Control  controller 
Kvc=50.0 

Vref^764.0 

D=  (Kvc*(Vref-Vc)  +  L*pIL)/E 

3. )  Multi-loop  controller 
hl=0.0056 
hN=1.6105 
hV=0.0058 

delILd  =  (IL-VcD/R) 
delVc=  (Vc-VcD) 

Vex  =  INTEG(delVc,0.0) 

!delD=  -hV*delVc  -  hI*deULd  -hN*Vcx 
!D  =  BOUND(1.0e-5,L0,  VcD/E  +  DelD) 
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1.)  ACSL  PROGRAM  WITH  CONSTANT  OUTPUT  VOLTAGE 

! 

program 

INITIAL 

MAXTERVAL  iflaxt  =  1  .Oe-5  !  "maximum  integration  step  size" 

MINTERVAL  mint  =  l.Oe-9 

ClNTERVALcint  =  5.0e-6  ! "data  communication  interval" 

algorithm  ialg  =  5  !  "integration  algorithm" 

!"4-R.K.  2nd,  5-RX  4th" 

NSTEPS  nstp  =  1 

CONSTANT  tstop  =  1.0  !  "stop  point  for  integration" 

CONSTANT  FreqT  =  5.0e+3  !  "switch  frequency" 

IL=126.8 

TD=  1.0/FreqT  !  "switching  time  length" 

iminx=0.1 

_  •pD  I  "dumby  variable  for  schedule  use" 

MINI  =  0.0  !  "value  ftw  comparison  purposes  only 

slope  =  1.0  !  "used  for  incrementing  WW  for  schedule" 

LOGICAL  TSCl,  TSOI  !  "TSCl  =  time  switch  closed  for  converter" 
TSCl  =  .true.  !  "TSOI  =  time  switch  open  for  converter" 

TSOI  =  .false. 

LOGICAL  CM,  NCM  !  "CM  =  continuious  mode" 

CM  =  .true.  !  "NCM  =  discontinuious  mode" 

NCM  =  .false. 
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LOGICAL  TD20,  TD2A  !  "TD20  =  time  D2  is  off  during  discontinuous  mode" 
TD20  =  .false.  !  "TD2A  =  time  D2  is  active  during  discontinuous  mode" 
TD2A  =  .true.  !  "D2  is  the  time  the  iduction  curent  becomes  zero" 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

CONSTANT  E=  850.0  !  "voltage  source  for  Buck" 

CONSTANT  Vc=  750.0  !  "voltage  across  the  load  of  Buck 

!  (note:  this  is  for  constant  load  R)" 

CONSTANT  L  =  1.35e-3  !  "Buck  inductance" 

CONSTANT  C  =  2600.0e-6  !"  buck  capacitance" 

CONSTANT  R  =  5.625  !  "load  resistance  of  buck" 

Imin  =  BOUND(0.0, 1.0e+34minx)  !  "current  limit  across  inductor 
inever  goes  below  zero" 

SCHEDULE  finalCMl  .xn.  IL  -  mini  !  "when  inductcx*  current  (TT.) 

!  "goes  to  zero  the  conveter  goes  discontinuious  " 

PROCEDURAL(D2  J),Imax4minx  =  Vc,E,L,R,TD)  !  "find  D  =  time  ratio  of 

iswitch  on/off' 

IF  (CM)  THEN  !  "continuous  mode" 

D  =  Vc/E 

Imax=  D*E*((1.0/R)+((1.0-D)*TD/(2*L))) 

Iminx=  D*E*((1.0/R)-((1.0-D)*TD/(2*L))) 

ELSE  !  "Discontinuous  mode" 

k  =  (Vc/E)**2.0/(2.0*(1.0  -  Vc/E)) 

D  =  sqrt((4.0*L*k)/(R*TD)) 

D2  =  (D/2.0)*(1.0  +  sqrt(1.0  +  (2.0/k))) 
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Imax  =  (Vc/L)*(D2-D)*TD 
Iminx  =  (Vc/R)  -  (VC/(2  *L))*(1.-D)*TD 
ENDIF 
END 

WT  =  INTEG(slope,  0.0) 

SCHEDULE  finallB  .xz.  WT  -  D*WW  !  "needed  to  initiate  TSCl  to  TSOI  at  D" 
SCHEDULE  final2B  .xz.  WT  -  WW  !  "needed  to  initiate  TSOI  to  TSCl 

!  at  new  cycle" 

SCHEDULE  finalD2  .xz.  WT  -  D2*WW  !  "needed  to  initiate  location  IL 
!  goes  to  zero  during  cycle" 

PROCEDURAUpIL  =  E,Vc,L,TSCl.TS01,CHNCM,Iniax,  diff) 

!  "finds  slope  of  IL" 

IF  ((TSCl)  AND.  (CM))  THEN 
pIL=  (E-Vc)/L 

ELSEIF  ((TSOI)  .AND.  (CM))  THEN 
pIL  =  -Vc/L 

EISEIF  ((TSCl)  AND.  (NCM))  THEN 
pIL=  (E  -  Vc)/L 

EISEIF  ((TSOI)  .AND.  (TD2A))  THEN 
pIL  =  -(Vc/L) 

ELSE 

pIL=0.0 

ENDIF 

END 

IL=  INTEG(pIL,0.01)  !  "finds  IL" 
power=  (Vc**2.0)/R 
END  !  "of  derivative" 

DISCRETE  finallB 
TSCl  =  (.false.) 

TSOI  =  (.true.) 
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END 

DISCRETE  final2B 
TSCl  =  Ctrae.) 
TSOI  =  (.false.) 
TD2A  =  (.true.) 
TD20  =  (.false.) 

WT  =  0.0 
IL=Imin 
END 

DISCRETE  finalCMl 
CM  =  (.false.) 

NCM  =  (.true.) 

END 

DISCRETE  finalCM2 
CM  =  (.true.) 

NCM  =  (.false.) 
TSCl  =  (.true.) 

TSOI  =  (.false.) 
TD2A  =  (.true.) 

TD20  =  (.false.) 
WT=0.0 
END 

DISCRETE  fmalD2 
TD20  =  (.true.) 
TD2A  =  (.false.) 
END 

END  !  "of  dynamic" 
END  !  "f  program" 


APPENDIX  D.  PSSCM  REDUCED-ORDER  MODEL  COPWLG^TION  WITH 
MUn-LOOP  CONTROLLER  WITH  FIXED  SOURCE 


PROGRAM 

INITIAL 

MAXTERVAL  maxt  =  l.Oe-5  ! "maximum  integration  step  size" 
MINTERVAL  mint  =  l.Oe-9 

QNTERVAL  cint  =  5.0e-6  !"data  communication  interval" 
ALGORITHM  ialg  =  5  !  "integration  algorithm" 

!"4~R.K.  2nd,  5~R.K.  4th" 


NSTEPS  nstp=l 

CONSTANT  tstop  =  1.0  !  "stop  point  fw  integration" 

CONSTANT  FreqT  =  5.0e+3  !  "switch  frequency" 

TD=  1.0/FreqT  !  "switching  time  length" 

WW  =  TD  !  "dumby  variable  Iot  schedule  use" 

slope  =  1.0  !  "used  fw  incrementing  WW  for  schedule" 

LOGICAL  TSCl,  TSOI  !  "TSCl  =  time  switch  closed  for  converter" 


TSCl  =  .true. 
TSOI  =  .false. 
Vcy=0.0 
Vc  =  750.0 
DelD=0.0 


!  'TSOI  =  time  switch  open  for  converter" 

!  "  i.c.  for  D  part  of  PID  controUw 
! "  initial  voltage  across  tire  load  of  Buck" 

! "  iiutial  perturbation  of  D  " 


END  !  "of  initial" 
DYNAMIC 


TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

CONSTANT  E=  850.0  !  "voltage  source  for  Buck" 

CONSTANT  Vcd=  750.0  !  "desired  voltage  across  the  load  of  Buck" 
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CONSTANT  L=  1.35e-3  !  "Buck  inductance" 

CONSTANT  C  =  2600.0e-6  !"  buck  capacitance" 

CONSTANT  R  =  5.625  !  "load  resistance  of  buck" 

WT  =  INTEG(slope,  0.0) 

SCHEDULE  finallB  .xz.  WT  -  D*WW  !  "needed  to  initiate  TSCl  to  TSOI  at  D" 
SCHEDULE  final2B  .xz.  WT  -  WW  !  "needed  to  initiate  TSOI  to  TSCl 

!  at  new  cycle" 

PROCEDURAL(pIL  =  E,Vc4.,TSCl,TS01)  !  "finds  slope  of  IL" 

IF  (TSCl)  THEN 
pIL=  (E-Vc)/L 
ELSEIF  (TSOI)  THEN 
pIL  =  -Vc/L 
ENDIF 
END 

IL=  BOUND(0.0,1000.0,  LIMINT(pIL,0.0,0.0, 1000.0))  !  "finds  IL" 

ir=  Vc/R 

ic=IL-ir 

Vc=  (l/c)*INTEG(ic,  0.0) 

!  Multiloop  controller 
hl=0.0056 
hN=1.6105 
hV=0.0058 
dellLd  =  (IL-VcD/R) 
delVc=  (Vc-VcD) 

Vex  =  INTEG(delVc,0.0) 

delD=  -hV*delVc-IiI*delILd -hN*Vcx 

D  =  BOUND(1.0e-5,1.0,  VcD/E  +  DelD) 

END  !  "of  derivative" 


87 


DISCRETE  finallB 
TSCl  =  (.false.) 
TSOI  =  (.true.) 
WTX  =  0.0 
END 

DISCRETE  final2B 
TSCl  =  (.true.) 
TSOI  =  (.false.) 
WT  =  0.0 
END 

END  !  "of  dynamic" 
end  !  "f  program" 


APPENDIX  E.  PSSCM  CONFIGURATION  WITH  MUTI-LOOP  CONTROLLER 

FOR  ACTIVE  SOURCE 


!"  This  program  runs  a  460v  50hp  induction  motor 

INCLUDE  'MACROS^mSarb.mac'  ! "  induction  machine  macro" 

PROGRAM 

INITIAL 

MAXTERVAL  maxt  =  1  .Oe-5  !  "maximum  integration  step  size" 

MINTERVAL  mint  =  l.Oe-7 

CENTER VAL  cint  =  4.0e-3  !  "data  communication  interval" 
ALGORITHM  ialg  =  5  !  "integration  algorithm" 

!"4-R.K.  2nd,  5-R.K.  4th" 

NSTEPS  nstp  =  1 

CONSTANT  tstop  =  3.0  !"stop  point  for  integration" 

CONSTANT  tpio3  =  2.094395 1  !  "two  pi  over  three" 
sqrt3  =  SQRT(3.0) 

CONSTANT  wrmpk  =  157.0  !"in  rad/sec  =  1500  RPM" 

CONSTANT  tacc  =  1.0  !"time  to  accelerate" 

CONSTANT  tcbinv  =  0.1 !  "time  at  which  inverter  (x>nnected" 

LOGICAL  cbinv!"(X)nnects  to  inverter/motor" 
cbinv  =  .FALSE.  !  "initially  discx>nnected” 

CONSTANT  FreqT  =  5.0e+3  !  "switch  frequency" 

CONSTANT  Vmag  =  890.0  !  "votage  magnitude  of  genterated  source' 

CONSTANT  L=  1.35e-3  !  SSCM  component  values 

CONSTANT  C  =  2600e-6 
CONSTANT  R  =  5.625 

spin  =  0.0  !  angle  for  rectifier  turn  on/off  dicxle 
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iS  =  0.0  !  ”  ic  current  through  converter  switch" 

irl  =  0.0  !"  ic  current  through  R" 

IL  =  0.0  ! "  ic  through  inductor" 

iavginv  =  0.0  !  "initial  average  inveter  current" 

VEic=  Vmag*(.95493)  ! "  ic  for  voltage  source  E" 

delD  =  0.0  ! "  initial  pertubation  of  duty  cycle" 

TD=  1.0/FreqT  !"time  duration  for  one  duty  cycle" 

=  TD  !"  variables  for  control  of  duty  cycle 

slope  =  1.0 

CONSTANT  VcD  =  750.0  !  "  inverter  input  voltage" 

!  "gives  max  volt  of  460" 

! "  rms  Une-to-line" 

D  =  .8824  !  "D=Vc/E" 

LOGICAL  TSCl,  TSOI  !  "TSCl  =  time  switch  closed  for  converter" 
TSCl  =  .true. !  "TSOI  =  time  switch  open  for  converter" 

TSOI  =  .false. 

LOGICAL  Multiloop 
Multiloop  =  .true. 

LOGICAL  Tlr,T2r,T3r,T4r,T5r,T6r  !"Recifier  logic" 

Tlr  =.true. 

T2r  =.true. 

T3r  =.false. 

T4r  =.false. 

T5r  =.false. 

T6r  =.false. 

!  "set  the  initial  ramp  up  rate  for  the  desired" 

!  "induct  motor  speed" 

IF  (cbinv)  THEN 
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pin  =  wrmpk/tacc 
ELSE 
pin  =  0.0 
ENDIF 

!  "initialize  inverter  switch  states  to  all  lower  closed" 

SA  =  .false. 

SB  =  .false. 

SC  =  .false. 

CONSTANT  TLrated  =  165.0  !  "rated  motor  load  torque" 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

CONSTANT  B1  =  0.0  ! "  friction  damping  coefficient" 

CONSTANT  J1  =  0.07  !  "  inertia  of  load  and  motor" 

!  "  assumed  equal  to  0.035  kg-m2" 

CONSTANT  wrmicl  =  0.0  !  "  rotor  speed  initial  cond." 
CONSTANT  thrmicl  =  0.0  !  "  rotor  position  i.c." 

CONSTANT  hyst  =  1.0  !"hysteresis  level  for  cxur  control" 
CONSTANT  idsestar  =  27.56  !"sets  sidre  =  360.48" 
CONSTANT  vOsl  =0.0 
CONSTANT  PI  =  3.141593 

CONSTANT  KKK  =  2.094395  !  "(2*pi)/3,  constant  for  3  phase" 
CONSTANT  Wb=  377.0 

CONST  ANT  Rlf=  0.236  !  "Filter  component  values 
CONSTANT  Xcf=  0.1332 
CONSTANT  Xlf=  4.725 

CONSTANT  kp  =  5.0! "proportional  gain  for  spd  control" 
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CONSTANT  ki  =  200.0!  "integral  gain  for  speed  control" 

!  "generates  the  desired  rotor  speed" 
wrmdes  =  INTEG(pin,  0.0) 

I  "schedule  when  ramp  up  of  speed  is  over" 

SCHEDULE  rampoff  .XP.  t-tacc-tcbinv 
1  "schedule  when  ramp  up  of  speed  begins" 

SCHEDULE  invon  .XP.  t-tcbinv 

Vab  =  Vmag*cos(spin  +  (KKK/4.0))  1  "3  phase  source" 

Vbc  =  Vmag*cos(spin  -  3*(KKK/4.0)) 

Vca  =  Vmag*cos(spin  +  5*(KKK/4.0)) 

CONSTANT  al  =  3600  I  "coefficinents  for  characteristic  equation" 
CONSTANT  b2  =  1.93e6 
CONSTANT  c3  =  3.9e8 

spin=INTEG(377.0, 0.0)  !"reg;ulates  the  rectifier  switching  speed" 

!  "produces  the  rectified  ouqjut  from  the  3  phase  source 
PROCEDURAL(Vr  =  Tlr,T2r,T3r,T4r,T5r,T6r,Vca,Vbc,Vab) 

F  ( (Tlr)  .AND.  (T2r))  THEN  !"Vr  is  rippled  DC  voltage" 
Vr=-Vca 

ELSEF  ((T2r)  .AND.  (T3r))  THEN 
Vr=Vbc 

ELSEF  ((T3r)  .AND.  (T4r))  THEN 
Vr=-Vab 

ELSEF  ((T4r)  .AND.  (T5r))  THEN 
Vr=Vca 

ELSEF  ((T5r)  .AND.  (T6r))  THEN 
Vr=-Vbc 

ELSEF  ( (Tlr)  .AND.  (T6r))  THEN 
Vr=  Vab 


92 


ENDIF 

END 

I’Turas  on  and  off  thyristors  for  rectifier 
SCHEDULE  fmallr  .XP.  spin  -  (2.0*PI  ) 

SCHEDULE  fmal2r  .XP.  spin  -  (PI/S.O ) 

SCHEDULE  finalSr  .XP.  spin  -  ((2.0/3.0)*PI  ) 

SCHEDULE  final4r  .XP.  spin  -  (PI ) 

SCHEDULE  fmalSr  .XP.  spin  -  ((4.0/3.0)*PI  ) 

SCHEDULE  fmal6r  .XP.  spin  -  ((5.0/3.0)*PI ) 

!  "LOW  PASS  FILTER" 

Lf=Xlf/Wb 

pir=(Vr-Vi-Rlf*irl)/Lf 
irl=LIMINT(pir.  0.0, 0.0,8000.0) 
pVi=(irl-iS)*(wb*Xcf) 

Vi=INTEG(pVi,  VEic)  !  "Smooth  dc  voltage" 

!  "BUCK  CONVERTER" 

E  =  Vi  !  "Input  voltage  to  Buck  converter" 

WT  =  INTEG(sIope,  0.0) !  "Ramp  control  for  control  frequency" 
SCHEDULE  fmallB  .xz.  WT  -  D*WW  {"Controls  switching  cycle" 
SCHEDULE  fmal2B  .xz.  WT  -  WW  !  "Controls  time  cycle" 
{"Control  to  determine  current  through  the  switch  and  inductor" 
PROCEDURAL(pIL,iS  =  E.Vc,L,TSCl,TS01,IL) 

IF(TSC1)  THEN 
pIL=  (E-Vc)/L 
iS  =  IL 

ELSEIF  (TSOI)  THEN 

pIL  =  -Vc/L 
iS  =  0.0 
ELSE 
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pIL=0.0 
iS  =  0.0 
ENDIF 
END 

IL=  BOUND(0.0, 1000.0,  LIMINT(pIL, 0.0, 0.0, 1000.0))  ["Inductor  current 

!  "must  not  go  below  zero" 

ir=Vc/R  ["current  through  the  risistor 

ic=  IL-iinv  -ir  ["Current  through  the  capacitor 

[ "  Vc  is  the  ou^ut  voltage  of  converter" 

Vc=  (l/c)*INTEG(ic,  1.95)  [  "750/(c=2.6e-3)  =  1.95  " 

Vcx  =  INTEG(delVc,0.0) 

Rim=  BOUND(1.0e-2,1000.0  ,Vc/(il  +  l.Oe-6)) 

[  Multiloop  controller 
El=  850.0 

D  =  BOUND(1.0e-5,1.0,  .8824  +  DelD) 

[  hI=(al*(L*R*C)  -  L)/(E1*  R  *C) 
hi=  .00571154 
[  hN=c3*(L*C)/El 
hN=  1.61047 

[  hV=(b2*(L*R*C)  -El*hl  -  R)/(E1*R) 
hV=  .00673618 

PROCEDURAL(delD,deIIld,delVc=Vcx,VcDJldL,VC) 

1 "  determines  change  in  current" 

IF  (multiloop)  THEN  1  "and  voltage  while  converter  is  on" 

delILd=  (IL  -  VcD/R  -  iavginv)  [ 
delVc  =  (Vc- VcD) 

delD  =  -hV*delVc  -  hI*delILd  -  hN*Vcx 
else 
delD  =0.0 
delILd  =  0.0 
delVc  =  0.0 
ENDIF 
END 

avgil=  INTEG(iinv,0.0) 

PROCEDlJRAL(Rinv,Rt=iavginv,cbinv,R)  [  "determines  resistance  with  I.M." 

IF  (abs(iavginv)  .LT.  .0001)  THEN 

Rinv=0.0 

Rt=R 

ELSEIF  (cbinv)  THEN 

Rinv=Vc/abs(iavginv) 

Rt  =  (Rinv*R)/(Rinv  +  R) 

else 

Rinv=0.0 
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Rt=R 

ENDBF 

END 

CONSTANT  wbl=  377.0 
!  " — invoke  induction  motor  MACRO" 
im3arb(  1  ,vqss  1 ,  vdss  1  ,vOs  1  ,wnn  1  ,wb  1 ,0.0,iqss  1  ddss  1  ,i0s  1 ,  & 
iqrsl,idrsUOrl,Tel."rsl  =0.087",  "irl  =0.228",  & 

"Xml  =  13.08",  "Xssl  =  13.382",  "Xrrl  =  13.382",  & 
"polesl=4.0") 

!  " — given  the  stationary  reference  frame  currents" 

!  "  establish  the  machine  currents" 
iasl  =iqssl 

ibsl  =-0.5*(iqssl+sqrt3*idssl) 
icsl  =-iasl  -  ibsl 

!  "determine  the  desired  slip  angle:  the  -  thr" 
thsldes  =  INTEG(wsldes,  0.0) 

!  "determine  the  desired  synchronous  ref.  frame  angle" 
thedes  =  thsldes  +  0.5*polesl*thrml 
iasdes  =  iqsestar*COS  (thedes)  +  idsestar*SIN(thedes) 
ibsdes  =  iqsestar*COS(thedes-q)io3)  +  idsestar*SIN(thedes-tpio3) 
icsdes  =  -iasdes  -  ibsdes 
! " — derive  a  simple  speed  control" 
speederr  =  wrmdes  -  wrml 
pxl  =  speederr 
xl  =  INTEG(pxl,  0.0) 
iqsestar  =  kp*speederr  +  ki*  xl 
!  " — develop  die  inverter  hysteresis  controls" 
PROCEDURAL(vasl,vbsl,vcsl,iinv=SA,SB,SC,iasl,ibsl,icsl,  & 
iasdes,ibsdes,icsdes,Vc,cbinv) 

!  "initialize  pole  voltages  to  zero" 
vap  =  0.0 
vbp  =  0.0 
vqj  =  0.0 

!  "initialize  currents  into  top  leg  switches  zero" 
iina  =  0.0 
iinb  =  0.0 
iinc  =  0.0 

!  "if  the  top  switch  is  closed  then  the  pole  voltage" 

!  "equals  the  input  voltage  and  the  current  through  " 

!  "that  switch  equals  the  phase  current" 
if  (SA)  then 
vap  =  Vc 
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iina  =  iasl 
endif 

if  (SB)  then 
vbp  =  Vc 
iinb  =  ibsl 
endif 

if  (SC)  then 
vq)  =  Vc 
iinc  =  icsl 
endif 

!  "calculate  the  current  into  the  inverter" 

IF  (cbinv)  THEN 

iinv  =  iina  +  iinb  +  iinc 

ELSE 

iinv  =  0.0 

ENDIF 

il=dinv 

!  "calculate  the  neutral  to  neg  input  line  voltage" 

vnp  =  (vap+vbp+vq))/3.0 

!  "calculate  the  phase  voltages" 

vasl  =  vap  -  vnp 

vbsl=vbp-vnp 

vcsl  =  vcp  -  vnp 

END 

I " —  mechanical  interface " 

IF  (cbinv)  THEN 

I " _ Given  the  machine  voltages,  determine  the  stat.  ref 

!  "frame  voltages" 
vqssl  =  vasl 

vdssl  =  (-vasl-2.0*vbsl)/sqrt3 
pwrml  =  (Tel  -  TL1)/J1 
pthrml  =  wrml 

! " — establish  the  reference  machine  currents" 

!  "calculate  the  desired  slip  frequency" 
wsldes  =  wbl*rrl*iqsestar/(Xrrl*idsestar) 

ELSE 
vqss  1  =  0.0 
vdssl  =  0.0 
pwrml  =  0.0 
pthrml  =  0.0 
wsldes  =  0.0 
ENDIF 

nwrml  =  wrml/188.5 
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pctTLl  =  0.11799  -  0.53915*abs(nwrml)  +  1.7*nwrml*nwrml  -  & 
0.28554*nwrm  1  *nwrm  1  *abs  (nwrm  1 ) 
PROCEDURAL(TLl=pctTLl  ,TLrated,nwrm  1) 

IF  (nwrml  .LT.  0.0)  THEN 
TLl  =  -pctTLl  *TLrated 
ELSE 

TLl  =  pctTLl  *TLrated 

ENDIF 

END 

wrml  =  INTEG(pwrml,  wrmicl) 
thrml  =  INTEG(pthrml,  thrmicl) 
iinvint  =  INTEG(iinv,  0.0) 

Pt=(Vc**2)/Rt 
Pinv  =Vc*iinv 
END  !  "of  derivative" 

!"  turns  on/off  the  rectifier  thyristors" 

DISCRETE  finallr 
T2r=  (.true.) 

T6r=(.false.) 
spin=  spin-2.0*PI 
end 

DISCRETE  final2r 
T3r=(.true.) 

Tlr=(.false.) 

end 

DISCRETE  fmaI3r 
T4r=(.true.) 

T2r=(.false.) 

end 

DISCRETE  final4r 
T5r  =(.true.) 

T3r=(.false.) 

end 

DISCRETE  fmalSr 
T6r  =(.true.) 

T4r  =(.false.) 
end 

DISCRETE  final6r 
Tlr  =(.true.) 

T5r=(.false.) 

end 

DISCRETE  final  IB!  "opens  converter  switch" 

TSCl  =  (.false.) 
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TSOI  =  (.true.) 
iS  =  0.0 
END 

DISCRETE  final2B  !"  closes  converter  switch 

TSCl  =  (.true.) 

TSOI  =  (.false.) 

WT  =  0.0 
END 

DISCRETE  average!"sets  up  measurement  for  avg  inveter" 
INTERVAL  tsample=  0.002 ["current" 
iavginv=  avgilAsample 
avgil=  0.0 
END 

DISCRETE  invon 

pin  =  wrmpk/tacc  !  "set  ramp  up  rate  for  des  spd" 
cbinv  =  .TRUE.  ["turn  inverter  on" 

END 

DISCRETE  rampoff 
pin  =  0.0 
wrmdes  =  157.0 
END 

DISCRETE  aveiinv 
INTERVAL  tsamp  =  0.002 
aveiin  =  iinvint/tsamp 
IF  (aveiin  .ne.  0.0)  THEN 
avezin  =  Vc^aveiin 
ENDIF 
iinvint  =  0.0 
END 

DISCRETE  swchng 
INTERVAL  tsamp2  =  0.0002 
!  "calculate  phase  current  errors" 
iaserr  =  iasdes  -  iasl 
ibserr  =  ibsdes  -  ibsl 
icserr  =  icsdes  -  icsl 

!  "if  phase  a  error  is  greater  than  +hyst  close  top  switch" 
if  (iaserr  .GT.  hyst)  then 
S  A  =  .true, 
endif 

!  "if  phase  a  error  is  less  than  -hyst  close  bot  switch" 
if  (iaserr  .LT.  -hyst)  then 
S  A  =  .false, 
endif 


98 


!  "if  neither  condition  is  true,  switch  status  stays  same" 
if  (ibserr  .GT.  hyst)  then 
SB  =  .true, 
endif 

if  (ibserr  .LT.  -hyst)  then 

SB  =  .false. 

endif 

if  (icserr  .GT.  hyst)  then 

SC  =  .true. 

endif 

if  (icserr  .LT.  -hyst)  then 

SC  =  . false. 

endif 

END 

END  !  "of  dynamic" 

END  !  "of  program" 
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APPENDIX  F.  MACRO  FOR  INDUCTION  MACHINE 


"z-induction  machine  identifier" 

ttii 

"  INPUTS" 

"vqs-q-axis  stator  voltage  in  the  warb  frame(V) 

"vds-d-axis  stator  voltage  in  the  warb  frame(V) 

"vOs-zero  seq.  stator  voltage  in  the  warb  frame(V) 

"wrm-rotor  mechanical  speed  (rad/sec) 

"wb-base  electrical  angular  velocity  (rad/sec)" 

"warb-speed  of  the  reference  frame  (rad/sec)" 

till 

"  OUTPUTS" 

"iqs-q-axis  stator  current  in  the  warb  frame  (A) 

"ids-d-axis  stator  current  in  the  warb  frame  (A) 

"iOs-zero  seq.  stator  current  in  the  warb  frame  (A) " 

"iqr-q-axis  rot-ref  current  in  the  warb  frame  (A) 

"idr-d-axis  rot-ref  current  in  the  warb  frame  (A) 

"iOr-zero  seq.  rot-ref  current  in  the  warb  frame  (A) " 

"Te-electromagnetic  torque,  positive  for  motor" 

"action  (N-m)" 

nil 

"  PARAMETERS" 

"rs&z-stator  winding  resistance  (ohms)" 

"rr&z-rotor-referred  winding  resistance  (ohms)" 

"Xm&z-stator  magnetizing  reactance  (ohms)" 

"Xss&z-stator  self  reactance  (ohms)" 

"sidsic&z  -d-axis  stator  flux  linkage/sec  initial  cond.  (V)" 

"siOsic&z  -0-seq.  stator  flux  linkage/sec  initial  cond.  (V)" 

"siqric&z  -q-axis  rot-ref  flux  linkage/sec  i.c.  (V) 

"sidric&z  -d-axis  rot-ref  flux  linkage/sec  i.c.  (V)" 

"siOric«&z  -0-seq.  rot-ref  flux  linkage/sec  i.c.  (V)" 

"D&z-constant  helpful  in  determining  currents  (ohms2) 

"Xls&z-stator  leakage  reactance  (ohms)" 

"Xlr&z-rotor-referred  leakage  reactance  (ohms)" 

"wr&z-rotor  electrical  speed  (rad/sec)" 

nn 

******************************************************************** 

MACRO  im3arb(z,vqs,vds,v0s,wrm,wb,warb,iqs,ids,i0s,iqr,idr,i0r,Te,  & 

prs,prr,pXm,pXss,pXrr  ,ppoles) 

INITIAL 
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CONSTANT  prs 
CONSTANT  prr 
CONSTANT  pXm 
CONSTANT  pXss 
CONSTANT  pXrr 
CONSTANT  ppoles 

!  "A  convenient  constant  for  detennining  the  currents" 

D&z  =  Xss&z*XiT&z  -  Xm&z*Xm&z 
!  "Establish  the  stator  and  rotor-refeired  leakage  reactance" 
Xls&z  =  Xss&z  -  Xm&z 
Xlr&z  =  Xrr&z  -  Xm&z 

!  "assign  the  flux  linkage  per  second  initial  conditions" 

CONSTANT  siqsic&z  =  0.0 

CONSTANT  sidsic&z  =  0.0 

CONSTANT  siOsic&z  =  0.0 

CONSTANT  siqric&z  =  0.0 

CONSTANT  sidric&z  =  0.0 

CONSTANT  siOric&z  =  0.0 

END 

!  "Compute  the  rotor  electrical  speed  given  the  mechanical  speed" 
wr&z  =  0.5*poles&z*wrm 

!  "Determine  the  machine  currents  from  the  state  varaiables" 

iqs  =  (XiT&z*siqs&z  -  Xm&z*siqr&z)/D&z 

ids  =  (Xrr&z*sids&z  -  Xm&z*sidr&z)/D&z 

iqr  =  (Xss&z*siqr&z  -  Xm&z*siqs&z)/D&z 

idr  =  (Xss&z*sidr&z  -  Xm&z*sids&z)/D&z 

!  "The  zero  sequence  currents" 
iOs  =  siOs&z/Xls&z 
iOr  =  siOr&z/Xlr&z 

!  "Calculate  the  derivative  of  the  flux  linkage  per  second" 
psiqs&z  =  -rs&z*wb*iqs  -  warb*sids&z  +  wb*vqs 
psids&z  =  -rs&z*wb*ids  +  warb*siqs&z  +  wb*vds 
psiqr&z  =  -rr&z*wb*iqr  -  (warb-wr&z)*sidr&z 
psidr&z  =  -rr&z*wb*idr  +  (warb-wr&z)*siqr&z 
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psiOs&z  =  -rs&z*wb*iOs  +  wb*vOs 
psiOr&z  =  -rr&z*wb*iOr 

!  "Integrate  the  derivatives  to  arrive  at  the  state  variable 

siqs&z  =  INTEG(psiqs&z,  siqsic&z) 

sids&z  =  INTEG(psids&z,  sidsic&z) 

siOs&z  =  INTEG(psiOs&z,  siOsic&z) 

siqr&z  =  INTEG(psiqr&z,  siqric&z) 

sidr&z  =  INTEGCpsidr&z,  sidric&z) 

siOr&z  =  INTEG(psiOr«&z,  siOric&z) 

!  "Compute  the  developed  electromagnetic  torque" 

Te  =  O.75*poles&z*(siqr&z*idr  -  sidr&z*iqr)/wb 

MACRO  END 


"Xrr&z-rotor-referred  self  reactance  (ohms)" 

"poles&z  -number  of  machine  poles 

ftM 

"INTERNAL  (STATE  OR  STATE  RELATED)" 
"siqs&z-q-axis  stator  flux  linkage/sec  warb  frame  (V) 
"sids&z-d-axis  stator  flux  linkage/sec  warb  frame(V)  " 
"siOs&z-O-seq.  stator  flux  linkage/sec  warb  frame(V) 
"siqr&z  -q-axis  rot-ref  flux  link/sec  warb  frame  (V)  " 

"sidr&z  -d-axis  rot-ref  flux  link/sec  warb  frame  (V) 
"siOr&z  -0-seq.  rot-ref  flux  link/sec  warb  frame  (V) 
"psiqs&z  -derivative  of  siqs&z  (V/sec)" 

"psids&z  -derivative  of  sids&z  (V/sec)" 

"psiOs&z  -derivative  of  siOs&z  (V/sec)" 

"psiqr&z  -derivative  of  siqr&z  (V/sec)" 

"psidr&z  -derivative  of  sidr&z  (V/sec) 

"psiOr&z  -derivative  of  siOr&z  (V/sec)" 

IMf 

"  internal  (NOT  STATE  RELATEDr 
"siqsic&z  -q-axis  stator  flux  linkage/sec  initial  cond.  (V)" 


APPENDIX  G.  STEAM  TURBINE  DRIVEN  SYNCHRONOUS  MACHINE  WITH 
EXCITOR  CONTOL-ACTIVE  SOURCE  INCLUDING  RECTIFIER  AND  FILTER 

1.)  ACSL  Program 

INCLUDE  ’MACROS/sm31kq.mac’ 

INCLUDE  'MACROS/Lncomconv3.mac' 

PROGRAM 

INITIAL 

MAXTER VAL  maxt  =  1  .Oe-7  !  "maximxim  integration  step  size" 

MINTERVAL  mint  =  l.Oe-9 

"To  accelerate  the  simulations  for  Rlarge  >  100" 

"start  with  ialg  =  5,  then  at  0.1  switch  to  Gears" 

"ialg  =2  with  maxt=1.0e-4" 

CINTERVAL  cint  =  1  .Oe-3  !  "data  communication  interval" 

ALGORITHM  ialg  =  5  !  "integration  algorithm" 

!"4-R.K.  2nd,  5-R.K.  4th" 

NSTEPS  nsti)  =  1 

CONSTANT  tstop  =  1.0  !"stop  point  for  integration" 

PARAMETER  (twopi  =  6.283 185307) 
al20  =  twopi/3.0 

"INITIAL  section  for  501  ssgt  model" 

CONSTANT  nref  =1.0  !  "per  unit  speed  ref  for  cntrl" 

CONSTANT  kcl  =  22.5  !"PI  controUer  gain" 

CONSTANT  tcl  =  0.55  !"PI  controUer  time  constant" 

CONSTANT  tfv  =  0.01,  tft  =  0.05 

CONSTANT  cgtl  =  1.3523,  cgt2  =  0.251,  wfos  =  0.23 

CONSTANT  cgn  =  0.5 

CONSTANT  qldic  =  0.0  !"Initial  prime  movr  torque" 

qloadi  =  qldic 

wfic  =  qloadi/cgtl  +  cgt2 

wfinic  =  wfic 

ucbic  =  (wfinic  -  wfos) 

"INi'ilAL  section  for  exciter  model" 

CONSTANT  keel  =  1.0,  teel  =  0.1 
CONSTANT  vrelul  =  8.4,  vrelU  =  0.0 
CONSTANT  kael  =  400.0,  tael  =  0.01,  kfel  =  0.01 
CONSTANT  talfel  =  0.15,  ta2fel  =  0.06,  taSfel  =  0.0 
CONSTANT  exfdic  =  1.0,  vbusic  =1.0 
vreli  =  exfdic/keel 
vfelxi  =  vreli 
vfelyi  =  kfel*vfelxi 
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vgref  =  vbusic  +  vreli/kael 

"INTIIAL  section  for  synchronous  machine" 

•• _ initial  flux  linkages  at  Ipu  voltage  and  zero  current" 

vgen  =1.0 

nwrmic2  =  1.0  !  "normalized  initial  speed" 

sikqric2  =  0.0 
sikdric2  =  exfdic 

sifdric2  =  exfdic  +  (0.136829/1. 768)*exfdic 
siqsric2  =  0.0 
sidsric2  =  1.0 
si0sric2  =  0.0 

"INmAL  section  for  the  lOOkW  buck  converter  load 
CONSTANT  VcD  =  550.0  ! "  Desired  output  voltage" 

CONSTANT  Rinl  =  7500.0 

CONSTANT  Cinl  =  3900.0e-6!  "input  filter  capacitance" 

CONSTANT  Linl  =  200.0e-6  !"input  filter  inductance" 

CONSTANT  ilinicl  =  0.0!"Linl  ic  current" 

CONSTANT  Vcapinil  =  0.0!"Cinl  ic  voltage" 

CONSTANT  dutymax  =  0.6!"max  duty  cycle" 

CONSTANT  Routl  =  5.625! "rated  load  resistance" 

CONSTANT  L  =  1 .35e-3  !  "Buck  inductance" 

CONSTANT  C  =  2600.0e-6  ! "  buck  capacitance" 

CONSTANT  Vbase  =  540.0!"base  voltage  3.125MVA  SM" 

"peak  phase  voltage" 

CONSTANT  Ibase  =  3.858e3!"base  current  3.125MVA  SM" 

"peak  phase  current" 

Vc  =0.0 

CONSTANT  tacc  =  0.5!  "time  to  bring  voltage  up" 

CONSTANT  tbuckon  =  2.0! "time  buck  1  turned  on" 

LOGICAL  buckup,Ibuckup, multiloop 
buckup  =  .True. 

Ibuckup  =.false. 
multiloop  =  .TRUE, 
ramprate  =  dutymax/tacc 
pduty  =  0.0 

LOGICAL  cbl,contmode 

cbl  =  .FALSE.  !"cb  connecting  buckl  to  bus" 

contmode  =  .TRUE. 

CONSTANT  FreqT  =  5.0e+3  !  "switch  frequency" 

TD=  1.0/FreqT  !  "switching  time  length" 

—  -pD  j  "dumby  variable  for  schedule  use" 

slope  =  1.0  !  "used  for  incrementing  WW  for  schedule" 

CONST ANT  k=25 .0  ! "  factor  for  estimateing  pole  location  for  M.L.  control' 
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LOGICAL  TSCl,  TSOI  !  "TSCl  =  time  switch  closed  for  converter' 
TSCl  =  .true.  !  "TSOI  =  time  switch  open  for  converter" 

TSOI  =  .false. 

END  !  "of  initial" 

DYNAMIC 

TERMT  (t  .GE.  (tstop-0.5*cint)) 

DERIVATIVE 

CONSTANT  wbl=  377.0 

" — determine  the  rotor  angle  for  machine  2" 

CONSTANT  deUc2  =  0.0 
pdel2  =  0.0 

del2  =  INTEG(pdel2.  deUc2) 

" — convert  the  synchronous  ref  frame  volt,  to  rotor  frame" 
vqsr2  =  COS(del2)*vqse2  -  SIN(del2)*vdse2 
vdsr2  =  SIN(del2)*vqse2  +  COS(del2)*vdse2 
CONSTANT  wb2  =  377.0 
sm3  Ikq(2,vqsr2,vdsr2,v0s2,vfdr2,wrm2,wb2,  & 
iqsr2,ikqr2,idsr2dfdr2,ikdr240s2,Te2,  & 

"rs2=0.00515","rkq2=0.0613",  & 

"rfd2=0.001 1 1  ","rkd2=0.023968","Xmq2=1.0",  & 
"Xq2=1.08","Xkq2=1.329787",  & 
"Xmd2=1.768","Xd2=1.848","Xkd2=2.101829",& 
"Xfd2=1.904829","poles2=2.0") 

" — convert  the  rotor  ref  frame  currents  to  e-frame" 
iqse2  =  COS(del2)*iqsr2  +  SIN(del2)*idsr2 
idse2  =  -SIN(del2)*iqsr2  +  COS(del2)*idsr2 
" — establish  the  rotor  angle  and  the  synchronous  angle" 
thrm2  =  INTEG(wrm2, 0.0) 

"find  e-frame  voltages  for  input  to  SM  macro" 
vqse2  =  2.0*(COS(thrm2)*vas2+COS(thrm2-al20)*vbs2  +  & 
COS(thrm2+al20)*vcs2)/3.0 

vdse2  =  2.0*(SIN(thnn2)*vas2+SIN(thrm2-al20)*vbs2  +  & 
SIN(thrm2+al20)*vcs2)/3.0 
v0s2  =  (vas2+vbs2+vcs2)/3.0 
" — convert  the  e-frame  currents  to  abc-variables" 

"inputs  to  Lncomconv  " 

ias2  =  COS(thrm2)*iqse2  +  SIN(thrm2)*idse2 

ibs2  =  COS(thrm2-al20)*iqse2  +  SIN(thrm2-al20)*idse2 

ics2  =  -ias2  -  ibs2 

"synchronous  machine  rotor  dynamics" 

CONSTANT  Hsm  =  2.137 

pnwrm2  =  (Te2-Tpm)/(2.0*Hsm) !  "norm,  deriv  of  speed" 
nwrm2  =  INTEG^nwrm2,  nwrmic2) 
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wnn2  =  wb2*nwrm2 

dngt  =  (Te2-Tpm)/(2.0*Hsm)  !"nonn.  deriv  of  speed" 
ngt  =  nwnn2 

•• — invoke  the  rectifier  model" 

Lncomconv(l, (-ias2),(-ibs2),(-ics2),Idcl,. true.,. true.,  & 
.true.,.true.,.true.,.true.,vdcl,vas2,vbs2,  & 
vcs2,"rlol=0.0001","rhil=1000.0  ) 

” — (Jc  link  dynamics" 

CONSTANT  rdcl  =  0.005 

CONSTANT  Ldcl  =  0.013 

pidcl  =  -rdcl*Idcl/Ldcl  +  (vdcl-vcap)/Ldcl 

Idcl  =  INTEG(pidcl,  0.0) 

CONSTANT  Ccap  =  0.1 
CONSTANT  rcap=10.0 

PROCEDlJRAL(pvcap=Idcl,vcap,iswl,Ibase,Ccap,rcap) 

JF  (cbl)  THEN 

"buck  convCTterl  connected" 

pvcap  =  (Idcl-vcap/rcap-iswl/Ibase)/Ccap 

else 

"buck  converter  1  disconnected" 
pvcap  =  (Idcl-vcapMap)/Ccap 
ENDIF 

END !  "of  procedural" 
vcap  =  INTEG(pvcap,  0.0) 

'• _ Calculate  the  generator  terminal  voltage  magnitude" 

vgen  =  SQRT(vqse2*vqse2  +  vdse2*vdse2) 

" — Implement  the  exciter  dynamic  model" 

eglxer  =  vgref  -  vgen  -  vfelb 

dvrel  =  (eglxer*kael  -  vrel)/tael 

vrel  =  LIMINT(dvrel,  vreli,  vrelU,  vrelul) 

vfelx  =  INTEG((vrel-vfelx)/ta2fel,  vfelxi) 

vfela  =  (ta3fel/ta2fel)*(vrel-vfelx)  +  vfelx 

vfelb  =  (kfel*vfela  -  vfely)/talfel 

vfely  =  INTEG(vfelb,  vfelyi) 

dexfd  =  (vrel  -  exfd*keel)/teel 

exfd  =  INTEG(dexfd,  exfdic) 

vfdr2  =  rfd2*exfd/Xmd2 

" — SSGT  Model" 

"input:  ngt  (per  unit  gen  speed)" 

"output:  Tpm  (power  turb  tork  pu  on  3125  kVA  base 

"Controller" 

nerr  =  nref  -  ngt  !  "normalized  speed  error" 

ductr  =  -kcl*dngt  +  (kcl/tcl)*nerr 
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uctr  =  INTEG(ductr,  uctric) 

"Fuel  Valve  and  Combustor  eqns" 
wfin  =  uctr  +  wfos 
REALPL(wfv  =  tfv,wfm,wfinic) 

REALPL(wf  =  tft,  wfvjwfic) 

qgtpu  =  cgtl*(wf-cgt2)4cgn*(1.0-ngt)  !"pu  on  2675  kVA" 
Tpm  =  -qgtpu*2675.0/3125.0 
" — invoke  the  rectifier  model" 
Lncomconv(l,(-ias2),(-ibs2),(-ics2),Idcl,. true., .true.,  & 
.true.,.true.,.true.,.true.,vdcl,vas2,vbs2,  & 
vcs2,"rlol=0.0001","rhil=1000.0") 

" — dc  linit  dynamics" 

CONSTANT  rdcl  =  0.005 
CONSTANT  Ldcl  =  0.013 
pidcl  =  -rdcl*Idcl/Ldcl  +  (vdcl-vcap)/Ldcl 
Idcl  =  INTEG(pidcl,  0.0) 

CONSTANT  Ccap  =  0.1 
CONSTANT  rcap=10.0 

PROCEDURAL(pvcap=Idcl ,  vcap,iswl  ,Ibase,Ccap,rcap) 

F  (cbl)  THEN 

"buck  converterl  connected" 

pvcap  =  (Idcl-vcapAcap-iswl/Ibase)/Ccap 

ELSE 

"buck  converterl  disconnected" 
pvcap  =  (Idcl-vcap/rcap)/Ccap 
ENDF 

END!  "of  procedural" 
vcap  =  INTEG(pvcap,  0.0) 

" — Calculate  the  generator  terminal  voltage  magnitude" 
vgen  =  SQRT(vqse2*vqse2  +  vdse2*vdse2) . 

" — Lnplement  the  exciter  dynamic  model" 

eglxer  =  vgref  -  vgen  -  vfelb 

dvrel  =  (eglxer*kael  -  vrel)/tael 

vrel  =  LIMINT(dvrel,  vreli,  vrelU,  vrelul) 

vfelx  =  INTEG((vrel-vfelx)/ta2fel,  vfelxi) 

vfela  =  (ta3fel/ta2fel)*(vrel-vfelx)  +  vfelx 

vfelb  =  Ocfel*vfela  -  vfely)/talfel 

vfely  =  INTEG(vfelb,  vfelyi) 

dexfd  =  (vrel  -  exfd*keel)/teel 

exfd  =  INTEG(dexfd,  exfdic) 

vfdr2  =  rfd2*exfd/Xmd2 

" — SSGT  Model" 

"input:  ngt  (per  unit  gen  speed)" 
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"output;  Tpm  (power  turb  tork  pu  on  3125  kVA  base" 
"Controller" 

nerr  =  nref  -  ngt  !  "normalized  speed  error" 
ductr  =  -kcl*dngt  +  (kcl/tcl)*nerr 
uctr  =  INTEG(ductr,  uctric) 

"Fuel  Valve  and  Combustor  eqns" 
wfin  =  uctr  +  wfos 
REALPL(wfv  =  tfv,wfm,wfinic) 

REALPL(wf  =  tft,  wfv.wfic) 

qgtpu  =  cgtl*(wf-cgt2)4cgn*(1.0-ngt)  !"pu  on  2675  kVA" 
Tpm  =  -qgtpu*2675.0/3125.0 
"implement  BUCK  converter  load  eqns" 

" — input  filter  eqns" 
vbuckin  =  Vbase*vcap 
pilinl  =  (vbuckin-vcapinl)/Linl 
iLinl  =INTEG(piLinl,  iLinicl) 

PROCEDURAUpvcapin  1  =iLin  1  .vcapin  1  ,isw  1  ,Cin  1  ,Rin  1  ,cb  1 ) 
JF  (cbl)  THEN 

pvcapinl  =  (iLinl-vcapinl/Rinl-iswl)/Cinl 

ELSE 

pvcapinl  =  (iLinl)/Cinl 
"pvcapinl  =  (iLinl-vcapinl/Rinl)/Cinl 
ENDIF 

END!  "of  procedural" 

vcapinl  =  INTEG(pvct^inl,  vcapinil) 

ioutl  =  Vc/Routl 

SCHEDULE  fullbuck  .XP.  t-tbuckon-tacc 
PROCEDURAL(pduty=cbl,buckup) 

IF  ((cbl)  .AND.  (buckup))  THEN 

pduty  =  ramprate 

ELSE 


pduty  =  0.0 

ENDIF 

END 


duty  =  INTEG(pduty,  0.0)  rrcz-ki 

SCHEDULE  fmallB  .xz.  WT  -  D*WW  I  "needed  to  initiate  TSCl  to  TSOI 

SCHEDULE  fmal2B  .xz.  WT  -  WW 
!  "needed  to  initiate  TSOI  to  TSCl  at  new  cycle" 

!  Multiloop  controller 
R=  Routl 

E=  Vcapinl  !  "voltage  source  for  Buck" 

El  =  875.0 


atD" 
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al  =3600.0 
b2  =  1.93e6 
c3=3.9e8 

M=  (al*(L*R*C)  -  L)/(E1*  R  *C) 
hN=  c3*(L*C)/El 

hV=  (b2*(L*R*C)  -El*hl  -  R)/(E1*R) 

PROCEDURAL(demd,delVc  =  VcD,Routl  JL,VC) 

!  "  determines  change  in  current" 

IF  ((cbl).and.  (Ibuckup))  THEN  ! "  voltage  while  converter  is  on" 
delILd  =  (IL-VcD/Routl) 
delVc=  (Vc-VcD) 
else 

delILd  =  0.0 
delVc  =  0.0 
ENDIF 
END 

Vex  =  INTEG(delVc,0.0) 

!"more  buck  equations" 

PROCEDURAL(DD  =  E,VcD,delD,Duty)  !  "  determines  d" 

IF  (Buckup)  THEN 
DD=Duty 
delD=0.0 

ELSEIF  ((Ibuckup)  .and.  (Multiloop))  THEN 
delD=  -hV*delVc  -  hI*delILd  -  hN*Vcx 
DD=VcD/E  +  delD 
ELSE 
delD=0.0 
DD=  0.6 
ENDIF 
END 

D  =  BOUND(1.0e-5,1.0,  DD) 

WT  =  INTEG(slope,0.0) 

PROCEDURAL(pIL,iswl  =  E,Vc,L,TSCl,TS01,IL,cbl  Jbase) 

!  "finds  slope  of  IL" 

JF  ((TSCl)  .and.  (cbl))  THEN 
pIL=  (E-Vc)/L 
iswl=  nvibase 

ELSEIF  ((TSOI)  .and.  (cbl))  THEN 
pIL  =  -Vc/L 
iswl=  0.0 
else 

pIL  =0.0 
isw  1=0.0 
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ENDIF 

END 

1L=  BOUND(0.0,1000.0,  LIMINT(pIL, 0.0,0.0,1000.0))  !  "finds  IL" 
ic=  EL-  (Vc^outl) 

Vc=  (l/c)*INTEG(ic,0.0) 

END  !  "of  derivative" 

" — schedule  end  of  up  ramp  for  the  duty  cycle 
DISCRETE  fuUbuck 
buckup  =  .false. 

Ibuckup  =.true. 
pduty  =  0.0 
duty  =  dutymax 
END!"offullbuck" 

DISCRETE  finallB 
TSCl  =  (.false.) 

TSOI  =  (.true.) 

WTX  =  0.0 
END 

DISCRETE  final2B 
TSCl  =  (.true.) 

TSOI  =  (.false.) 

WT  =  0.0 
END 

END  !  "of  dynamic" 

©ID  !  "of  program" 

?***”^TiI^********************************************************" 

M*****************************************************************" 

"Author:  John  G.  Ciezki" 

"Last  Revised:  26  Apr  94" 

"Developed  For:  Personal  Use" 

"  MACRO  Title:  LncomconvS 
"DESCRIPTION:" 

"ASSUMPTIONS:" 

"VARIABLE  DEFINITION" 

••  CONCATENATION" 

"z-3-phase  line-commutated  converter  identifier" 

"  INPUTS" 

"ia-a-phase  current  referenced  positive  into" 

"the  converter  (A)" 

"ib-b-phase  current  referenced  positive  into" 

"the  converter  (A)" 
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"ic-c-phase  current  referenced  positive  into" 

"the  converter  (A)" 

"Idc-dc-side  current  referenced  positive  out  of 
"the  converter  (A)" 

"vgl-gate  signal  for  thyristor  1  ( T  or  F )" 

"vg2-gate  signal  for  thyristor  2  ( T  or  F )" 

"vg3-gate  signal  for  thyristor  3  ( T  or  F )" 

"vg4-gate  signal  for  thyristor  4  ( T  or  F )" 

"vg5-gate  signal  for  thyristor  5  ( T  or  F )" 

"vg6-gate  signal  for  thyristor  6  ( T  or  F )" 

"  OUTPUTS" 

"vdc-dc-side  voltage  positive  node  on  top  (V)" 

"va-a-phase  line-to-neutral  voltage  (V)" 

"vb-b-phase  line-to-neutral  voltage  (V)" 

"vc-c-phase  line-to-neutral  voltage  (V)" 

"  PARAMETERS" 

"rlo&z-the  value  of  the  bimodal  thyristor  resistance" 

"during  conduction  (ohms)" 

"rhi&z-the  value  of  the  bimodal  thyristor  resistance" 

"during  blocking  (ohms)  " 

"INTERNAL  (STATE  OR  STATE  RELATED)" 

"none" 

"  INTERNAL  (NOT  STATE  RELATED)" 

"j&z-index  for  going  through  all  possible  conv  states" 
"index&z  -records  index  of  proper  conduction  configuration" 
"cl&z-possible  thyristor  1  conduction  states  " 

"0  =  blocking,  1  =  conducting" 

"c2&z-possible  thyristor  2  conduction  states  " 

"0  =  blocking,  1  =  conducting" 
c2&z-possible  thyristor  2  conduction  states  " 

"0  =  blocking,  1  =  conducting" 

"c3&z-possible  thyristor  3  conduction  states  " 

"0  =  blocking,  1  =  conducting" 

"c4&z-possible  thyristor  4  conduction  states  " 

"0  =  blocking,  1  =  conducting" 

"c5&z-possible  thyristor  5  conduction  states  " 

"0  =  blocking,  1  =  conducting" 

"c6«&z-possible  thyristor  6  conduction  states  " 

"0  =  blocking,  1  =  conducting" 

"  cl&z  0  00  000000000000  etc" 

"  c2&z  000000000  0  0  0  0  0  0  etc" 

"  c3&z  000000001  1  1  1  1  1  1  etc" 

"  c4&z .0  00011110000  1  1  1  etc" 
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"  c5&z  0011001100  1  1  00  1  etc" 

••  c6&z  010101010101010  etc" 

"index  12  3  4  5  67  8  91011  12  131415 
"vl&z-logical  conduction  status  of  thyristor  1 " 

"F  =  not  conducting,  T  =  conducting" 

"v2&z-logical  conduction  status  of  thyristor  2" 

"F  =  not  conducting,  T  =  conducting" 

"v3«&z-logical  conduction  status  of  thyristor  3" 

"F  =  not  conducting,  T  =  conducting" 

"v4&z-logical  conduction  status  of  thyristor  4" 

"F  =  not  conducting,  T  =  conducting" 

"v5&z-logical  conduction  status  of  thyristor  5" 

T  =  not  conducting,  T  =  conducting" 

"v6&z-logical  conduction  status  of  thyristor  6" 

T  =  not  conducting,  T  =  conducting" 

"vallowl&z  -logical  status  of  whether  thyristor  1  can  be" 

"conducting  current,  F  =  no,  T  =  yes' 

"vallow2&z  -logical  status  of  whether  thyristor  2  can  be" 
"conducting  current,  F  =  no,  T  =  yes' 

"vallow3&z  -logical  status  of  whether  thyristor  3  can  be" 
"conducting  current,  F  =  no,  T  =  yes" 

"vallow4&z  -logical  status  of  whether  thyristor  4  can  be" 
"conducting  current,  F  =  no,  T  =  yes 
"vallow5&z  -logical  status  of  whether  thyristor  5  can  be" 
"conducting  current,  F  =  no,  T  =  yes" 

"vallow6&z  -logical  status  of  whether  thyristor  6  can  be" 
"conducting  current,  F  =  no,  T  =  yes' 

"VALVEl&z(64)  -matrix  of  the  possible  conduction  states  for 
"thyristor  1  for  the  64  possible  conv  configs." 

"VALVE2&z(64)  -matrix  of  the  possible  conduction  states  for" 

"thyristor  2  for  the  64  possible  conv  configs." 

"VALVE3&z(64)  -matrix  of  the  possible  conduction  states  for 
"thyristor  3  for  the  64  possible  conv  configs." 

"VALVE4&z(64)  -matrix  of  the  possible  conduction  states  for" 

"thyristor  4  for  the  64  possible  conv  configs." 

"VALVE5&z(64)  -matrix  of  the  possible  conduction  states  for 

"thyristor  6  for  the  64  possible  conv  configs." 

"CIDC&z(64)  -matrix  of  impedance  values  which  when  multiplied 

"by  Idc  yields  a  component  of  vdc  (ohms) 

CIA&z(64)  -matrix  of  impedance  values  which  when  multiplied 

"by  ia  yields  a  component  of  vdc  (ohms)" 

"C3B&z(64)  -matrix  of  impedance  values  which  when  multiplied 
"by  ib  yields  a  component  of  vdc  (ohms) 
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"vdc  =  CE)C&z()*Idc  +  CIA&z()*ia  +  CIB&zO*ib" 
"rl&z-the  bimodal  resistance  representing  thyristor  1 " 
"r2&z-the  bimodal  resistance  representing  thyristor  2  " 
"r3&z-the  bimodal  resistance  representing  thyristor  3  " 
"r4&z-the  bimodal  resistance  representing  thyristor  4 " 
*'r5&z-the  bimodal  resistance  representing  thyristor  5  " 
"r6&z-the  bimodal  resistance  representing  thyristor  6 " 
"all  in  (ohms)" 

"  Idc— >" 

•» _  ti 

"III  +" 

"III" 

"rl&zr3&zr5&z  " 

"  ia  ™>l  ib  — >1  ic  — >1" 

"III  vdc" 

"r4&zr6&zr2&z  " 

"III" 


" _  II 

II  gii 

"vdcmax&z  -initially  set  to  alrge  negative  number,  it  will  " 

"be  the  largest  possible  dc  voltage  for  the " 

"given  currents  ia,  ib,  ic  and  Idc  (V)" 

"vag&z-the  a-phase  to  ground  (-  terminal)  voltage  (V)" 

"vbg&z-the  b-phase  to  groimd  (-  terminal)  voltage  (V)" 

"vcg&z-the  c-phase  to  ground  (-  terminal)  voltage  (V)" 

"vng&z-the  neutral-to-ground  voltage  for  the  wye-conn." 

"ac  system  (V)" 

!"!!!!!!!!!!!!!!  MACRO  DEFINITION  BEGINS  HERE  !!!!!!!!!!!!!!!!!!!" 

MACRO  Lncomconv(z,ia,ib,ic,Idc,vg  1  ,vg2,vg3,vg4,vg5,vg6,vdc,va,vb,vc,  & 
prlo,  prhi) 

INITIAL 

CONSTANT  prlo 
CONSTANT  prhi 

"Define  indexing  integers" 

INTEGER  j&z,  index&z 

INTEGER  cl&z,  c2&z,  c3«S:z,  c4&z,  c5&z,  c6&z 


113 


"Define  conduction  logical  variables" 

LOGICAL  vl&z,  v2&z,  v3&z,  v4&z,  v5&z,  v6&z 
"Define  future  conduction  logical  variables" 

LOGICAL  vallowl&z,  vallow2&z,  vallow3&z 
LOGICAL  vallow4&z,  vallow5&z,  vallow6&z 
"Define  matrices  of  logic  states  for  64  possible  conv.  configs 
LOGICAL  VALVEl&z(64) 

LOGICAL  VALVE2&z(64) 

LOGICAL  VALVE3&z(64) 

LOGICAL  VALVE4&z(64) 

LOGICAL  VALVE5&z(64) 

LOGICAL  VALVE6&z(64) 

"Dimension  impedance  vectors  needed  to  represent  impedance  seen 
"  by  currents  ia,  ib  and  Idc  through  to  vdc 
DIMENSION  CIA&z(64) 

DIMENSION  CIB&z(64) 

DIMENSION  CIDC&z(64) 

"Loop  through  all  64  possible  conduction  states" 

DO  lO&z  cl&z  =  0,1 
DO  lO&z  c2&z  =  0,1 
DO  lO&z  c3&z  =  0,1 
DO  lO&z  c4&z  =  0,1 
DO  lO&z  c5&z  =  0,1 
DO  lO&z  c6&z  =  0,1 

"If  the  state  is  =  1  then  set  the  conduction  logic  T  " 

vl&z  =  cl&z  .EQ.  1 

v2&z  =  c2&z  .EQ.  1 

v3&z  =  c3&z  .EQ.  1 

v4&z  =  c4&z  .EQ.  1 

v5&z  =  c5&z  .EQ.  1 

v6&z  =  c6&z  .EQ.  1 

"If  conducting,  set  bimodal  resistance  to  rlo&z" 

"  otherwise  set  it  to  rhi&z" 
rl&z  =  RSW(vl&z,  rlo&z,  rhi&z) 
r2&z  =  RSW(v2&z,  rlo&z,  rhi&z) 
r3&z  =  RSW(v3&z,  rlo&z,  rhi&z) 
r4&z  =  RSW(v4&z,  rlo&z,  rhi&z) 
r5&z  =  RSW(v5&z,  rlo&z,  rhi&z) 
r6&z  =  RSW(v6&z,  rlo&z,  rhi&z) 

"Keep  running  index  of  possible  configs  ->  1  to  64" 
index&z  =  1  +  c6&z  +  2*(c5&z  +  2*(c4&z  +  2*(c3&z  +  & 
2*(c2&z  +  2*cl&z)))) 

"If  the  thryristor  is  conducting  at  a  particular  value" 
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"of  index,  then  set  the  vector  location  to  T" 
VALVEl&z(index&z)  =  vl&z 
VALVE2&z(index&z)  =  v2&z 
VALVE3&z(index&z)  =  v3&z 
VALVE4&z(index&z)  =  v4&z 
VALVES &z(index&z)  =  vS&z 
VALVE6&z(index&z)  =  v6&z 
"Compute  the  impedance  needed  to  compute  the  Idc" 
"contribution  to  the  dc  voltage  vdc" 


1  1  1" 


"rl+r4  r3  +  r6  r2  +  r5" 

CIDC&z(index&z)  =  -1.0/(1.0/(rl&z+r4&z)+1.0/(r3&z+r6&z)  & 
+  1.0/(r2&z+r5&z)) 

"Compute  the  impedance  needed  to  compute  the  ia  " 

"contribution  to  the  dc  voltage  vdc" 

CIA&z(index&z)  =  (r4&z/(r2&z+r4&z)-rl&2/(rl&z+r5&z))*  & 
((rl&z+r5&z)*(r2&z+r4&z)/  & 
(rl&z+r2&z+r4&z+r5&z))*  & 

(r3&z+r6&z)/  & 

(r3&z+r6&z+(r2&z+r5&z)*(rl«&z+r4&z)/  & 
(r2&z+r5&z+rl&z+r4&z)) 

"Compute  the  impedance  needed  to  compute  the  ib  " 

"contribution  to  the  dc  voltage  vdc" 

CIB&z(index&z)  =  (r6&z/(r2&z+r6&z)-r3&?/(r3&z+r5«fez))*  & 
((r3&z+r5&z)*(r2&z+r6&z)/  & 

(r3&z+r5&z+r2&z+r6&z))*  & 

(rl«Sfcz+r4&z)/  & 

(rl&z+r4&z+(r2&z+r5&z)*(r3&z+r6&z)/  & 
(r2&z+r5&z+r3&z+r6&z)) 
lO&z..  CONTINUE 

"Initialize  valve  conduction  to  no  valves  conducting" 

vl&z  =  .FALSE. 

v2&z  =  .FALSE. 

v3&z  =  .FALSE. 

v4&z  =  .FALSE. 

v5&z=. FALSE. 

v6&z=. FALSE. 

END!  "OF  INITIAL" 

"Establish  the  vdc  by  evaluating  it  for  each  of  the  64" 
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"converter  configurations  with  ia,  ib,  ic,  and  Idc  and  using 
"the  maximum  value  found" 

PROCEDURAI/vdc,  rl&z,  r2&z,  r3&z,  r4&z,  rS&z,  r6&z  -  & 
vl&z,  v2&z,  v3&z,  v4&z,  v5&z,  v6&z,  & 
vgl,  vg2,  vg3,  vg4,  vg5,  vg6,  & 

ia,  ib,  Idc)  ^ 

"A  switch  may  conduct  if  it  has  been  conducting  or 

"  if  a  gate  signal  is  applied" 

vallowl&z  =  vl&z  .OR.  vgl 

vallow2&z  =  v2&z  .OR.  vg2 

vallow3&z  =  v3&z  .OR.  vg3 

vallow4&z  =  v4&z  .OR.  vg4 

vallow5&z  =  v5&z  .OR.  vg5 

vallow6&z  =  v6&z  .OR.  vg6 

"very  small  value  to  begin  iteration  of  comparisons" 

vdcmax&z  =  -1.0e30 

index&z  =  0 

"cycle  through  the  64  pssible  configurations" 

DO  20&Z  j&z  =  1,64 

"Check  if  the  assumed  configuration  is  a  " 

"conduction  configuration  or  whether  the  valve" 

"could  be  clamped  on" 

IF  (((.NOT.  VALVEl&z&O&z))  .OR.  vallowl&z )  .AND.  & 
((.NOT.  VALVE2&z&(j&z))  .OR.  vallow2&z )  .AND.  & 
((.NOT.  VALVE3&z&(j&z))  .OR.  vaUow3&z )  .AND.  & 
((.NOT.  VALVE4&z&(j&z))  .OR.  vallow4&z )  .AND.  & 
((.NOT.  VALVE5&z&(j&z))  .OR.  vallow5&z )  .AND.  & 
((.NOT.  VALVE6&z&(j&z))  .OR.  vallow6&z ))  THEN 
"If  the  configuration  is  valid,  then" 

"compute  the  dc  voltage" 

vdc  =  CIA&z(j&z)*ia  +  CIB&zO&z)*ib  +  & 
CIDC&z(j&z)*Idc 

"compare  the  new  value  of  the  dc  voltage  with" 

"the  previous  maximum" 

JP  (vdc  .GT.  vdcmax&z)  THEN 

"If  the  new  value  is  greater  set  it" 

"equal  to  the  maximum,  and  store  the" 

"index  indicating  that  configuration" 
vdcmax&z  =  vdc 
index&z  =  j&z 

END  IF 
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END  IF 

20&Z..  CONTINUE 

"the  dc  voltage  equals  the  maximum  voltage  foimd  among" 

"the  valid  conduction  configurations" 
vdc  =  vdcmax&z 

"set  the  valve  conduction  thyristor  according  to  the" 

"desired  configuration  index" 
vl&z  =  VALVEl&z&(index&z) 
v2&z  =  VALVE2&z&(index&z) 
v3&z  =  VALVE3&z&(index&z) 
v4&z  =  VALVE4&z&(index&z) 
v5&z  =  VALVES &z&(index&z) 
v6&z  =  VALVE6&z&(index&z) 

"set  the  actual  value  of  the  bimodal  resistors" 

rl&z  =  RSW(vl&z,  rlo&z,  rhi&z) 

r2&z  =  RSW(v2&z,  rlo&z,  rhi&z) 

r3&z  =  RSW(v3&z,  rlo&z,  rhi&z) 

r4&z  =  RSW(v4&z,  rlo&z,  rhi&z) 

rS&z  =  RSW(v5&z,  rlo&z,  rhi&z) 

r6&z  =  RSW(v6&z,  rlo&z,  rhi&z) 

END  !  "OF  PROCEDURAL" 

"given  the  desired  resistor  values  and  dc  voltage," 

"determine  the  phase  to  ground  voltages" 
vag&z  =  r4&z*(vdc+rl&z*ia)/(rl&z+r4&z) 
vbg&z  =  r6&z*(vdc+r3&z*ib)/(r3&z+r6&z) 
vcg&z  =  r2&z*(vdc+r5&z*ic)/(r2&z+r5&z) 

"establish  the  neutral-to-ground  voltage" 

vng&z  =  (vag&z  +  vbg&z  +  vcg&z)/3.0 

"calculate  the  phase  voltages" 

va  =  vag&z  -  vng&z 

vb  =  vbg&z  -  vng&z 

VC  =  vcg&z  -  vng&z 

MACRO  END 

3)  Sm31kq  Macro 

"♦************************************************************M?****" 

"Author:  John  G.  Ciezki" 

"Last  Revised:  26  Apr  94" 

"Developed  For:  Personal  Use" 

"  MACRO  Title:  sm31kq.mac  " 

"DESCRIPTION:" 

"ASSUMPTIONS:" 
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"VARIABLE  DEFINITION" 

"  CONCATENATION" 

"z-three  phase  synchronous  machine  identifier" 

"  INPUTS" 

"vqsr-q-axis  stator  voltage,  rotor  reference  fi-ame  (V)" 
"vdsr-d-axis  stator  voltage,  rotor  reference  fi-ame  (V)" 
"vOs-zero  sequence  stator  voltage  (V)" 

"vfdr-referred  field  voltage,  rot  ref  frame  (V)" 

"wrm-rotor  mechanical  speed  (rad/sec)" 

"wb-base  electrical  angular  velocity  (rad/sec)" 

"  OUTPUTS" 

"iqsr-q-axis  stator  current,  rot  ref  frame  (A)" 

"idsr-d-axis  stator  current,  rot  ref  firame  (A)" 

"ikqlr-#l  q-axis  referred  damper  current  (A)" 

"ikq2r-#2  q-axis  referred  damper  current  (A)" 

"ikdr-d-axis  referred  damper  current  (A)" 

"ifdr-d-axis  referred  field  current  (A)" 

"iOs-zero  sequence  stator  current  (A)" 

"Te-electromagnetic  torque,  positive  for" 

"motor  action  (N-m)" 

"  PARAMETERS" 

"rs&z-stator  phase  winding  resistance  (ohms)" 

"rkql&z-#l  q-axis  referred  damper  resistance  (ohms)" 
"rkq2&z-#2  q-axis  referred  damper  resistance  (ohms)" 
"rkd&z-d-axis  referred  damper  resistance  (ohms)" 
"rfd&z-d-axis  referred  field  resistance  (ohms) " 
"Xmq&z-q-axis  magnetizing  reactance  (ohms)" 
"Xq&z-q-axis  reactance  (ohms) 

"Xkql&z-#1  q-axis  referred  damper  self  inductance  (ohms)" 
"Xkq2&z-#2  q-axis  referred  damper  self  inductance  (ohms)" 
"Xmd&z-d-axis  magnetizing  reactance  (ohms)" 
"Xd&z-d-axis  reactance  (ohms)" 

"Xkd«fez-d-axis  referred  damper  self  inductance  (ohms)" 
"Xfd&z-d-axis  referred  field  self  inductance  (ohms)" 
"poles&z  -number  of  machine  poles" 

"INTERNAL  (STATE  OR  STATE  RELATED)" 
"siqsr&z-q-axis  stator  flux  linkage/sec,  rotor  reference " 
"fi-ame  (V)" 

"sidsr&z-d-axis  stator  flux  linkage/sec,  rotor  reference  " 
"fi-ame  (V)" 

"siOsr&z-zero  sequence  stator  flux  linkage/sec  (V) 
"sikqlr&z  -#1  q-axis  referred  damper  flux  linkage/sec, 
"rotor  reference  frame  (V)" 
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"sikq2r&z  -#2  q-axis  referred  damper  flux  linkage/sec,  " 

"rotor  reference  frame  (V)" 

"sikdr&z-d-axis  referred  damper  flux  linkage/sec,  rotor  " 

"reference  frame  (V)" 

"sifdr&z-d-axis  referred  field  flux  linkage/sec,  rotor  " 

"reference  frame  (V)" 

"psiqsr&z  -derivative  of  siqsr&z  (V/sec)" 

"psidsr&z  -derivative  of  sidsr&z  (V/sec)" 

"psiOsr&z  -derivative  of  siOsr&z  (V/sec)" 

"psikqlr&z  -derivative  of  sikqlr&z  (V/sec)" 

"psikq2r&z  -derivative  of  sikq2r&z  (V/sec)" 

"psikdr&z  -derivative  of  sikdr«&z  (V/sec)" 

"psifdr&z  -derivative  of  sifdr&z  (V/sec)" 

"  INTERNAL  (NOT  STATE  RELATED)" 

"Xd3&z-constant  helpful  in  determining  d-axis  currents  " 
"(ohms**3)" 

"adll&z-l  idsr  I  I  adll&z  adl2&z  adl3&z  I  Isidsr&zl  " 

"adl2&z-l  ifdr  I  =  I  adl2&z  ad22&z  ad23&z  I  Isifdr&zl  " 

"adl3&z-l  ikdr  I  I  adl3&z  ad23&z  ad33&z  I  Isikdr&zl  " 

"ad22&z-" 

"ad23&z-" 

"ad33&z-" 

"Xq3&z-constant  helpful  in  determining  q-axis  currents  " 
"(ohms**3)" 

"bql  l&z-l  iqsr  I  I  bql  l&z  bql2&z  bql3&z  I  lsiqsr«&z  I " 

"bql2&z-l  ikqlrl  =  I  bql2&z  bq22&z  bq23«&z  I  Isikqlr&zl  " 
"bql3&z-l  ikq2rl  I  bql3&z  bq23«&z  bq33«&z  I  Isikq2r&zl  " 
"bq22«fez-" 

"bq23&z-" 

"bq33&z-" 

"siqsric&z  -q-axis  flux  linkage/sec  initial  condition  (V)" 

"sidsric&z  -d-axis  flux  linkage/sec  initial  condition  (V)" 

"siOsric&z  -zero  seq.  flux  linkage/sec  initial  condition  (V)" 
"sikqlric&z  -#1  q-axis  referred  damper  flux  linkage/s  i.c.(V)" 
"sikq2ric&z  -#2  q-axis  referred  damper  flux  linkage/s  i.c.(V)" 
"sikdric&z  -d-axis  referred  damper  flux  linkage/sec  i.c.  (V)" 
"sifdric&z  -d-axis  referred  field  flux  linkage/sec  i.c.  (V)" 

"wr&z-rotor  electrical  angular  velocity  (rad/sec)" 

!"!!!!!!!!!!!!!!  MACRO  DEFINITION  BEGINS  HERE  !!!!!!!!!!!!!!!!!!!" 
MACRO  sm31kq(z,vqsr,vdsr,v0s,vfdr,wrm,wb,iqsr,ikqr,  & 
idsr, ifdr, ikdr, iOs.Te,  & 
prs,prkq,prfd,prkd,  & 
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pXmq,pXq,pXkq,pXmd,pXd,pXkd,pXfd, ,  & 

ppoles) 

INITIAL 

CONSTANT  prs 

CONSTANT  prkq 

CONSTANT  prfd 

CONSTANT  prkd 

CONSTANT  pXmq 

CONSTANT  pXq 

CONSTANT  pXkq 

CONSTANT  pXmd 

CONSTANT  pXd 

CONSTANT  pXkd 

CONSTANT  pXfd 

CONSTANT  ppoles 

"Determine  the  stator  leakage  reactance" 

Xls&z  =  Xq&z  -  Xmq&z 

"A  convenient  constant  for  determining  q-axis  currents 

Xq2&z  =  Xq&z*Xkq&z  -  Xmq&z*Xmq&z 

"determine  the  array  elements  for  relating  q-axis  flux 

"linkage  per  sec  vars  to  q-axis  currents" 

bql  l&z  =  Xkq&z/Xq2&z 

bql2&z  =  -Xmq&z/Xq2&z 

bq21&z  =  -Xmq&z/Xq2&z 

bq22&z  =  Xq&z/Xq2«Srz 

"A  convenient  constant  for  determining  d-axis  currents" 

Xd3&z  =  Xmd&z*Xmd&z*(2.0*Xmd&z-Xd&z-Xkd&z- 

Xfd&z)+Xd&z*Xfd&z*Xkd&z 

"determine  the  array  elements  for  relating  d-axis  flux  linkage" 

"per  second  variables  to  the  d-axis  currents" 
adll&z  =  (Xfd&z*Xkd&z  -  Xmd&z*Xmd&z)/Xd3&z 
adl2«&z  =  (Xmd&z*Xmd&z  -  Xkd&z*Xmd&z)/Xd3&z 
adl3&z  =  (Xmd&z*Xmd&z  -  Xfd&z*Xmd&z)/Xd3&z 
ad22&z  =  (Xd&z*Xkd&z  -  Xmd&z*Xmd&z)/Xd3&z 
ad23&z  =  (Xmd&z*Xmd&z  -  Xd&z*Xmd&z)/Xd3&z 
ad33&z  =  (Xd&z*Xfd&z  -  Xmd&z*Xmd&z)/Xd3&z 
"Set  the  state  variable  initial  conditions" 

" — set  in  main  program" 

"CONSTANT  siqsric&z  =  0.0" 

"CONSTANT  siiric&z  =0.0" 

"CONSTANT  siOsric&z  =  0.0" 

"CONSTANT  sikqric&z  =0.0" 

"CONSTANT  sifdricifez  =0.0" 
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"CONSTANT  sikdric&z  =0.0" 

END 

"Determine  the  rotor  electrical  speed  from  the  mechanical  speed" 
wr&z  =  0.5*poles&z*wrm 

"Determine  the  rotor  electrical  speed  from  the  mechanical  speed" 
wr&z  =  0.5*poles&z*wrm 

"Calculate  the  q-axis  currents  from  q-axis  flux  linkages/sec" 
iqsr  =  bql  l&z*siqsr&z  +  bql2&z*sikqr&z 
ikqr  =  bq21&z*siqsr&z  +  bq22&z*sikqr&z 
"Calculate  the  d-axis  currents  from  d-axis  flux  linkages/sec" 
idsr  =  adll&z*sidsr&z  +  adl2&z*sifdr&z  +  adl3&z*sikdr&z 
ifdr  =  adl2&z*sidsr&z  +  ad22&z*sifdr&z  +  ad23&z*sikdr&z 
ikdr  =  adl3&z*sidsr&z  +  ad23&z*sifdr&z  +  ad33&z*sikdr&z 
"Calculate  the  0  sequence  stator  current" 
iOs  =  siOsr&z/Xls&z 

"Establish  the  state  variable  derivative  equations" 

psiqsr&z  =  -rs&z*wb*iqsr  -  wr&z*sidsr&z  +  wb*vqsr 

psidsr&z  =  -rs&z*wb*idsr  +  wr&z*siqsr&z  +  wb*vdsr 

psiOsr&z  =  -rs&z*wb*iOs  +  wb*vOs 

psikqr&z  =  -rkq&z*wb*ikqr 

psifdr&z  =  -rfd&z*wb*ifdr  +  wb*vfdr 

psikdr&z  =  -rkd&z*wb*ikdr 

"Integrate  the  state  variables" 

siqsr&z  =  INTEG(psiqsr&z,  siqsric&z) 

sidsr&z  =  INTEG(psidsr&z,  sidsric&z) 

siOsr&z  =  INTEG(psiOsr&z,  siOsric&z) 

sikqr&z  =  INTEG(psikqr&z,  sikqric&z) 

sifdr&z  =  INTEG(psifdr&z,  sifdric&z) 

sikdr&z  =  INTEG(psikdr&z,  sikdric&z) 

"Compute  the  developed  electromagnetic  torque" 

"Te  =  0.75*poles&z*(sidsr&z*iqsr  -  siqsr&z*idsr)/wb" 

Te  =  sidsr&z*iqsr  -  siqsr&z*idsr  !"in  per  unit" 

MACRO  END 
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